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Abstract 

We consider the two logarithmic strain measures 

Wiso = ||dev n log U\\ = ||dev„ log VF T F\\ and cu vo i = |tr(log C/)| = |tr(log VF T F)\ = |log(det[/)| , 

which are isotropic invariants of the Hencky strain tensor log U, and show that they can be uniquely 
characterized by purely geometric methods based on the geodesic distance on the general linear group 
GL(n). Here, F is the deformation gradient, U = V F T F is the right Biot-stretch tensor, log denotes 
the principal matrix logarithm, || . || is the Frobenius matrix norm, tr is the trace operator and dev„ X = 
X — - tr(X) • 1 is the n-dimensional deviator of A' £ R riX ”. This characterization identifies the Hencky (or 
true) strain tensor as the natural nonlinear extension of the linear (infinitesimal) strain tensor e = sym Xu, 
which is the symmetric part of the displacement gradient Xu, and reveals a close geometric relation 
between the classical quadratic isotropic energy potential 

||dev„ sym Vu|| 2 + ^ [tr(sym Vtt)] 2 = ft ||dev n e|| 2 + ^ [tr(e)] 2 

in linear elasticity and the geometrically nonlinear quadratic isotropic Hencky energy 

M ||dev n log U || 2 + b [tr(log U)f = n w? so + ^ a; 2 ol , 

where fr is the shear modulus and k denotes the bulk modulus. Our deduction involves a new fundamental 
logarithmic minimization property of the orthogonal polar factor R, where F = R U is the polar decom¬ 
position of F. We also contrast our approach with prior attempts to establish the logarithmic Hencky 
strain tensor directly as the preferred strain tensor in nonlinear isotropic elasticity. 
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1 Introduction 

1.1 What’s in a strain? 

The concept of strain is of fundamental importance in elasticity theory. In linearized elasticity, one assumes 
that the Cauchy stress tensor a is a linear function of the symmetric infinitesimal strain tensor 

£ = symVu = sym(V</j — 1) = sym (F — 1), 

where </?: fl —> R" is the deformation of an elastic body with a given reference configuration fl C R n , 
x) = x + u(x ) with the displacement it, F = X7<p is the deformation gradient^ sym Vu = |(Vu + (Vw) T ) 
is the symmetric part of the displacement gradient Vw and 1 € GL + (n) is the identity tensor in the group 
of invertible tensors with positive determinant. In geometrically nonlinear elasticity models, it is no longer 
necessary to postulate a linear connection between some stress and some strain. However, nonlinear strain 

1 Although F is widely known as the deformation “gradient”, F = S/p = Dip actually denotes the first derivative (or the 
Jacobian matrix) of the deformation p. 
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tensors are often used in order to simplify the stress response function, and many constitutive l aws are 
expressed in terms of linear relations between certain strains and stresset0 15j 16, 24] (cf. Appendix A.2 for 
examples). 

There are many different definitions of what exactly the term “strain” encompasses: while Truesdell and 
Toupin [204] p. 268] consider “any uniquely invertible isotropic second order tensor function of \t he right 
Cauchy-Green deformation tensor C = F T F}” to be a strain tensor, it is commonly assumed 106 p. 230] 
(cf. 107. 108 23, 159]) that a (material or LagrangiaiQ strain takes the form of a primary matrix function 


of the right Biot-stretch tensor U = V F T F of the deformation gradient F £ GL + (n ), i.e. an isotropic tensor 
function E: Sym + ( n ) —> Sym(n) from the set of positive definite tensors to the set of symmetric tensors of 
the form 


E(U) = J2 c ( A Tg 


for 


^ = E A * 


(i.i) 


with a scale function e: (0,oo) —► R, where (g> denotes the tensor product, A i are the eigenvalues and ei are 
the corresponding eigenvectors of U. However, there is no consensus on the exact conditions for the scale 
function e; Hill (cf. EnG 459] and |l08l p. 14]) requires e to be “suitably smooth” and monotone with 
e(l) = 0 and e'(l) = 1, whereas Ogden 162, p. 118] also requires e to be infinitely differentiable and e' > 0 
to hold on all of (0, 00 ). 

The general idea underlying these definitions is clear: strain is a measure of deformation (i.e. the change 
in form and size) of a body with respect to a chosen (arbitrary) reference configuration. Furthermore, the 
strain of the deformation gradient F £ GL + (n) should correspond only to the non-rotational part of F. In 
particular, the strain must vanish if and only if A is a pure rotation, i.e. if and only if F £ SO(n), where 
SO(n) = {Q £ GL(n) | Q T Q = 1, detQ = 1} denotes the special orthogonal group. This ensures that the 
only strain-free deformations are rigid body movements: 


F t F = 1 


Vip(x) = F(x ) = R(x) £ SO(n) 
ip(x) = Qx + b for some fixed Q £ SO(n), b £ 


( 1 . 2 ) 


where the last implication is due to the rigidity 174 inequality ||Curli?|| 2 > c + ||Vi?|| 2 for R £ SO(n) (with 


a constant c + >0), cf. 151 . A similar connection between vanishing strain and rigid body movements holds 
for linear elasticity: if e = 0 for the linearized strain e = sym Vit, then u is an infinitesimal rigid displacement 
of the form 

u{x) = Ax + b with fixed A s so(n), b £ R n , 

where so (n) = {A £ R nx ” : A r = — A} denotes the space of skew symmetric matrices. This is due to the 


inequality ||Curl A|| 2 > c + ||VA|| 2 for A £ so(n), cf. 151 


In the following, we will use the term strain tensor (or, more precisely, material strain tensor ) to refer 
to an injective isotropic tensor function U 1 —> E(U) of the right Biot-stretch tensor U mapping Sym + (n) to 
Sym(n) with 


and 


E(Q t UQ) = Q t E(U)Q 
E(U) = 0 U = t: 


for all Q £ 0(n) 


(isotropy) 


where O(n) = {Q £ GL(n) | Q T Q = 1} is the orthogonal group and 1 denotes the identity tensor. In 
particular, these conditions ensure that 1 = U = 'JF T F if and only if F £ SO(n). Note that we do not 
require the mapping to be of the form (1.11. 


2 In a short note [32] , R. Brannon observes that “usually, a researcher will select the strain measure for which the stress-strain 
curve is most linear”. In the same spirit, Bruhns [33[ p. 147] states that u we should [.. .] always use the logarithmic Hencky 


strain measure in the description of finite deformations. Truesdell and Noll [203| p. 347] explain: “ Various authors [. . .] have 
suggested that we should select the strain [tensor] afresh for each material in order to get a simple form of constitutive equation. 
[.. .] Every invertible stress relation T = f(B) for an isotropic elastic material is linear, trivially, in an appropriately defined, 
particular strain [tensor f(B)]” 

3 Similarly, a spatial or Eulerian strain tensor E(V) depends on the left Biot-stretch tensor V = VFF T (cf. [ 74 ] ). 
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Among the most common examples of material strain tensors used in nonlinear elasticity is the Seth-Hill 
family^ 190 

1) :rel\{0} 

log U 


E r {U) = 


: r = 0 


(1.3) 


of material strain tensor^ which includes the Biot strain 
tensor Ey 2 (U ) = U — 1, the Green-Lagrangian strain ten¬ 
sor E^U) = \(C- 1) = \(U 2 - 1), where C = F T F = U 2 
is the right Cauchy-Green deformation tensor, the (ma¬ 
terial) Almansi strain tensor 0 E-i(U) = \{t -C” 1 ) 
and the (material) Hencky strain tensor Eq(U ) = log U, 
where log: Sym + (n) —> Syrn(n) is the principal matrix 
logarithm 105 p. 20] on the set Sym + (n) of positive def¬ 
inite symmetric matrices. The Hencky (or logarithmic) 
strain tensor has often been considered the natural or true 
strain in nonlinear elasticity 198 197, 75} 88]. It is also 
of great imp ortance t o so-called h ypoela stic models, as is 
(cf. Section 4.2.1) A very useful 


discussed in 210 76 



approximation of the material Hencky strain tensor was 
given by Bazant [17 (cf. (HHUHI): 

E y2 (U) := \ {E 1/2 {U) + E_ 1/2 {U)] = 

(1.4) 


Figure 1: Scale functions e r ,e r associated with the strain 
tensors E r and E r = ( ( E r — E— r ) via eigenvalue 
A. 


Additional motivations of the logarithmic strain tensor were also given by Vallee 205 206 , Rougee 182 


p. 302] and Murphy 142 . An extensive overview of the properties of the logarithmic strain tensor and its 

and 


applications can be found in 209 


154 


All strain tensors, by the definition employed here, can be seen as equivalent: since the mapping U K > E(U) 
is injective, for every pair E , E' of strain tensors there exists a mapping ip : Sym(n) —► Sym(n) such that 
E'(U) = ip(E(U)) for all U £ Sym + (n). Therefore, every constitutive law of elasticity can - in principle - be 
expressed in terms of any strain tensoiQand no strain tensor can be inherently superior to any other strain 
tensor]^] Note that this invertibility property also holds if the definition by Hill or Ogden is used: if the strain 

4 Note that log U = lim -)—[U 2r — 1). Many more examples of strain tensors used throughout history can be found in 47 

i — ^0 z r 

and [58], cf. [27] p. 132], 

5 The corresponding family of spatial strain tensors 


Er(V) = 


— J Jr (y 2r — 1 ) -.r^O 


logy 


: r = 0 


includes the Almansi-Hamel strain tensor ^1/2 (V) = V — 1 as well as the Euler-Almansi strain tensor E— i(V) = — B - 1 ), 

where B = FF T = V 2 is the Finger tensor 69 . 

6 Bruhns [37[ p. 41-42] emphasizes the advantages of the Hencky strain tensor over the other Seth-Hill strain tensors in the 
one-dimensional case: “The significant advantage of this logarithmic (Hencky) measure lies in the fact that it tends to infinity 
as F tends to zero, thus in a very natural way bounding the regime of applicability to the case F > 0. This behavior can also 
be observed for strain [tensors] with negative exponent n. Compared with the latter, however, the logarithmic measure also goes 
to infinity as F does, whereas it is evident that for negative values of n the strain ( F n — 1)] is bound to the limit 

All measures with positive values of n including the Green strain share the reasonable property of the logarithmic strain for 
F going to infinity. For F going to zero, however, these measures arrive at finite values for the specific strains, e.g. at —i 
for n = 2, which would mean that interpreted from physics a total compression of the rod (to zero length) is related to a finite 
value of the strain. This awkward result would not agree with our observation - at least what concerns the behavior of solid 
materials” 

7 According to Truesdell and Toupin [204| p. 268], .. any [tensor] sufficient to determine the directions of the principal axes 

of strain and the magnitude of the principal stretches may be employed and is fully general”. Truesdell and Noll [203[ p. 348] 
argue that there “is no basis in experiment or logic for supposing nature prefers one strain [tensor] to another”. 

8 Nevertheless, u [in] spite of this equivalence, one strain [tensor] may present definite technical advantages over another one” 
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is given via a scale function e, the strict monotonicity of e implies that the mapping U K > E(U) is strictly 
monotone 130 , i.e. 

(E(U 1 )-E(U 2 ),U 1 -U 2 )>0 

for all U i , U 2 G Sym + (n) with U\ 7 ^ U 2 , where (X,Y) = tr(X T y) denotes the Frobenius inner product on 
Sym(n) and tr(X) = y)" =1 is the trace of X G K raxrl . This monotonicity in turn ensures that the mapping 
U i->- E(U) is injective. 

In contrast to strain or strain tensor, we use the term strain measure to refer to a nonnegative real¬ 
valued function to: GL + (n) —> [0,oo) depending on the deformation gradient which vanishes if and only if F 
is a pure rotation, i.e. 00 (F) = 0 if and only if F G SO(n). 

Note that the terms “strain tensor” and “strain measure” are sometimes used interchangeably in the 
literature (e.g. 108 159| ). A simple example of a strain measure in the above sense is the mapping F 1-7 
||if(v / F T F)|| of F to an orthogonally invariant norm of any strain tensor E. 

There is a close connection between strain measures and energy functions in isotropic hyperelasticity: an 
isotropic energy potential [84] is a function W depending on the deformation gradient F such that 


W(F) >0, 
W(QF) = W(F) 
W(FQ) = W(F) 

for all F G GL + (n), Q G SO (n) and 


W(F) =0 if and only if F G SO(n). 


(normalization) 
(frame-indifference) 
(material symmetry: isotropy) 

(stress-free reference configuration) 


While every such energy function can be taken as a strain measure, many additional conditions for “proper” 


energy functions are discussed in the literature, such as constitutive inequalities 202 106 107 11 44 127 


generalized convexity conditions 10 13 or monotonicity conditions to ensure that “stress increases with 


strain” 154 Section 2.2]. Apart from that, the main difference between strain measures and energy functions 
is that the former are purely mathematical expressions used to quantitatively assess the extent of strain in 
a deformation, whereas the latter postulate some physical behaviour of materials in a condensed form: an 
elastic energy potential, interpreted as the elastic energy per unit volume in the undeformed configuration, 
induces a specific stress response functiorj^ and therefore completely determines the physical behaviour of 
the modelled hyperelastic material. The connection between “natural” strain measures and energy functions 
will be further discussed later on. 

In particular, we will be interested in energy potentials which can be expressed in terms of certain 
strain measures. Note carefully that, in contrast to strain tensors, strain measures cannot simply be used 
interchangeably: for two different strain measures (as defined above) 001 , 002 , there is generally no function 
/: R + —► R + such that 002 (F) = f(ooi(F)) for all F G GL + (n). Compared to “full” strain tensors, this can 
be interpreted as an unavoidable loss of information for strain measures (which are only scalar quantities). 

Sometimes a strain measure is employed only for a particular kind of deformation. For example, on the 
group of simple shear deformations (in a fixed plane) consisting of all F 1 G GL + (3) of the form 


F 7 = 


/I 7 0\ 
0 i 0 

Vo 0 


we could consider the mappings 


1 


F ^2~< ' 


171 


or 


7 € 


F 7 4 = In ( J 

7 73 l 2 




47| p. 467]. For example, there is one and only one spatial strain tensor E together with a unique objective and corotational 

d D _ d lo s • 


rate -Sp such that dD 

at 


dt E = sym (FF 1 ) = D. Here, ^ i s the logarithmic rate, D is the unique rate of stretching and E 


4.2.1 


and 


is the spatial Hencky strain tensor Eq = log V ; cf. Section 

9 The specific elasticity tensor further depends on the particular cJ 
constitutive law. 


a stress tensor in which to express the 
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We will come back to these partial strain measures in Section [3.2| 

In the following we consider the question of what strain measures are appropriate for the theory of 
nonlinear isotropic elasticity. Since, by our definition, a strain measure attains zero if and only if F £ SO(n), 
a simple geometric approach is to consider a distance function on the group GL + (n) of admissible deformation 
gradients, i.e. a function dist: GL + (n) x GL + (n) —► [0, oo) with dist(A,l?) = dist(R,A) which satisfies the 
triangle inequality and vanishes if and only if its arguments are identical^] Such a distance function induces 
a “natural” strain measure on GL + (n) by means of the distance to the special orthogonal group SO(n): 

w(F) := dist(F, SO(n)) := inf dist (F,Q). (1.5) 

QeSO(n) 

In this way, the search for an appropriate strain measure reduces to the task of finding a natural, intrinsic 
distance function on GL + (n). 


1.2 The search for appropriate strain measures 

The remainder of this article is dedicated to this task: after some simple (Euclidean) examples in Section [2j 
we consider the geodesic distance on GL + (n) in Section [3j Our main result is stated in Theorem 3.3 if the 
distance on GL + (n) is induced by a left-GL(n)-invariant, right-0(n)-invariant Riemannian metric on GL(ro), 
then the distance of F £ GL + (n) to SO(ra) is given by 

distg eod (E,SO(n)) = dist 2 geod (F,R) = p ||dev n log Uf + ^ [tr(logt/)] 2 , 

where F = RU with U = VF T F £ Sym + (?r) and R £ SO(n) is the polar decomposition of F. Section [i] also 
contains some additional remarks and corollaries which further expand upon this Riemannian strain measure. 

In Section |4j we discuss a number of different approaches towards motivating the use of logarithmic strain 
measures and strain tensors, whereas applications of our results and further research topics are indicated in 
Section 0 


Our main result (Theorem 3.3) has previously been announced in a Comptes Rendus Mecanique article 


147 as well as in Proceedings in Applied Mathematics and Mechanics 148 


The idea for this paper was conceived in late 2006. However, a number of technical difficulties had to be 
overcome (cf. [29) [l56j [U8j [l29j [145] ) in order to prove our results. The completion of this article might have 
taken more time than was originally foreseen, but we adhere to the old German saying: Gut Ding will Weile 
haben. 


2 Euclidean strain measures 


2.1 The Euclidean strain measure in linear isotropic elasticity 

An approach similar to the definition of strain measures via distance functions on GL + (n), as stated in 


equation (1.5), can be employed in linearized elasticity theory: let <p(x) = x + u(x) with the displacement u. 


Then the infinitesimal strain measure may be obtained by taking the distance of the displacement gradient 
Vu £ R nxrl to the set of linearized rotations so(n) = {A £ R nxn : A T = —A}, which is the vector spac^j 
of skew symmetric matrices. An obvious choice for a distance measure on the linear space R nxn = R” of 
n x n-matrices is the Euclidean distance induced by the canonical Frobenius norm 


I AT 


= Mxt X ) = ^'£ x * 


i,j =1 


10 A distance function is more commonly known as a metric of a metric space. The term “distance” is used here and throughout 
the article in order to avoid confusion with the Riemannian metric introduced later on. 

11 Note that so(n) also corresponds to the Lie algebra of the special orthogonal group SO(n). 
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We use the more general weighted norm defined by 


\\ X \\l^, K = A 4 l|dev„ symX|| 2 + p c ||skewX|| 2 + | [tr(X)] 2 , 


n, He, K > 0 , 


( 2 . 1 ) 


which separately weights the deviatoric (or trace free ) symmetric part dev n sym X = sym X — ^ tr(sym X) ■ 1, 
the spherical part. r tr(X) • 1, and the skew symmetric part skew X = |(X — X T ) of X; note that HXH^ K = 
||X|| for n = p c = 1, k = and that || • || jU)A4c K is induced by the inner product^ 2 ] 

(X, Y)n tlic>K = p (dev„ sym A', dev„ symfi) + p c (skew A, skew A) + § tr(X) tr(A) (2.2) 

on R nx ", where (A, V) = tr(X T A) denotes the canonical inner product. In fact, every isotropic inner product 
on R nx ”, i.e. every inner product (q-jiso with 


(Q t XQ, Q T YQ) iso = (X,Y)i. 


for all A,7 £ R nxn and all Q £ O (n), is of the form ( |2.2[ ), cf. 50 . The suggestive choice of variables p and 
k, which represent the shear modulus and the bulk modulus , respectively, will prove to be justified later on. 
The remaining parameter p c will be called the spin modulus. 

Of course, the element of best approximation in so(n ) to Vm with respect to the weighted Euclidean 
distance dist Euclid, Y) = \\X — Y || K is given by the associated orthogonal projection of Xu to 
so (n), cf. Figure [2j Since so(n) and the space Sym(n) of symmetric matrices are orthogonal with respect 
to (•, •) jU>/iC)K , this projection is given by the continuum rotation , i.e. the skew symmetric part skewVw = 
t(Vu — (Xu) T ) of Xu, the axial vector of which is curlw. Thus the distance iW* 


dist Euclid ^, MciK (Vu,so(n)) : = inf ||Vw-A|| 

AEso(n) h 


= || — skew V-u,|| = 11 sym Vu|| 


(2.3) 


We therefore find 


dist Euc iid !AllAlc , K (VM, so(n)) = ||sym Xu\\l^ K 

= p ||dev„ sym Vm|| 2 + ^ [tr(symVu)] 2 
= p ||dev„ e|| 2 + | [tr(e)] 2 = W Vln {Xu) 


for the linear strain tensor e = symVu, which is the quadratic isotropic elastic energy, i.e. the canonical 
model of isotropic linear elasticity with 


<t = -DvuWii n (Vu) = 2p clev„£ + Ktr(e) • 1. 


(2.4) 


12 The family ( [272] ) of inner products on R nXn is based on the Cartan-orthogonal decomposition 

0[(n) = ^ sl(n) PI Sym(n)^ ® so(n) ® R • 1 

of the Lie algebra Ql(n) = R nXn . Here, sl(n) = {X E gl(n) | trX = 0} denotes the Lie algebra corresponding to the special 
linear group SL(rz) = {A E GL(n) | det A = 1}. 

13 The distance can also be computed directly: since 

|| Viz — = fi 11dev n sym(Viz A)|| 2 +/z c 11skew(Viz - A)|| 2 + ^ [tr(Viz — A)] 2 

= /z ||dev n sym Viz|| 2 + /z c ||(skew Viz) — A|| 2 + — [tr(Viz)] 2 

for all A E so(n), the infimum inf || Viz — A|| K = (i ||dev n sym Viz|| 2 + ^ [tr(Viz)] 2 is obviously uniquely attained at 


A = skew Viz. 


Aeso(n) 
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Figure 2: The Euclidean distance dist|, uclid K (Vu,so(n)) = // |jdev n e|| 2 + [tr(e)] 2 of Vu to so(n) in R nXn in the infinites¬ 

imal strain setting. The strain tensor e = symVu is orthogonal to the infinitesimal continuum rotation skew Vu. 


This shows the aforementioned close connection of the energy potential to geometrically motivated measures 
of strain. Note also that the so computed distance to so(n) is independent of the parameter /z c , the spin 
modulus , weighting the skew-symmetric part in the quadratic form (2.1|. We will encounter the (lack of) 
influence of the parameter fi c subsequently again. 

Furthermore, this approach motivates the symmetric part e = sym Vu of the displacement gradient as 
the strain tensor in the linear case: instead of postulating that our strain measure should depend only on e, 
the above computations deductively characterize e as the infinitesimal strain tensor from simple geometric 
assumptions alone. 


2.2 The Euclidean strain measure in nonlinear isotropic elasticity 

In order to obtain a strain measure in the geometrically nonlinear case, we must compute the distance 

dist(V<^, SO(n)) = dist(F, SO(n)) = inf dist (F,Q) 

QeSO(n) 

of the deformation gradient F = V<p £ GL + (n) to the actual set of pure rotations SO(n) C GL + (n). It 
is therefore necessary to choose a distance function on GL + (n); an obvious choice is the restriction of the 
Euclidean distance on R" xrl to GL + (n). For the canonical Frobenius norm ||.||, the Euclidean distance 
between F, P £ GL + (n) is 


distEuclid (F, P ) = II F - P|| = yJtT[{F-P) T (F-P)]. 

Now let Q e SO(n). Since || . || is orthogonally invariant, i.e. ||QX|| = ||XQ|| = ||X|| for all X € R nxn , 
Q £ O(n), we find 

dist Euclid (E,Q) = \\F — Q\\ = \\Q t (F-Q)\\ = ||Q T E-1||. (2.5) 

Thus the computation of the strain measure induced by the Euclidean distance on GL + (n) reduces to the 
matrix nearness problem [104 

distEuclid (F, SO(ra)) = inf ||F - Q|| = min \\Q T F - t\\ . 

QESO(n) QeSO(n) 


By a well-known optimality result discovered by Giuseppe Grioli [82| (cf. [1501 H [m} [31]), also called 


“Grioli’s Theorem” by Truesdell and Toupin 204 p. 290], this minimum is attained for the orthogonal polar 
factor R. 

















Theorem 2.1 (Grioli’s Theorem 


82 


150 


204 


). Let F £ GL + (n). Then 


min \\Q t F - 1|| = || R t F - 1|| = \\VWf - 1|| = \\U - 
QeSO(n) 


where F = RU is the polar decomposition of F with R = polar(T') £ SO(n) and U = y/F T F £ Sym + (n). 
The minimum is uniquely attained at the orthogonal polar factor R. 

Remark 2.2. The minimization property stated in Theorem |2.1| is equivalent to [132] 

max tr (Q T F) = max ( Q T F , 1) = (R T F, 1) = (U. 1) . □ 

QeSO(n) QeSO(n) 

Thus for nonlinear elasticity, the restriction of the Euclidean distance to GL + ( n ) yields the strain measure 

dist Euclid (-F, SO(n)) = ||t/- 1||. 


In analogy to the linear case, we obtain 



Figure 3: The “flat” interpretation of Gib (n) C R” x n endowed with the Euclidean distance. Note that | F — R\\ = \\R(U— 1)|| = 
\\U — 1|| by orthogonal invariance of the Frobenius norm, where F = RU is the polar decomposition of F. 


distLciid(^SO(n)) = \\U-tf = \\E 1/2 \\ 2 , 


( 2 . 6 ) 


where Ey 2 = U — t is the Biot strain tensor. 
Venant-KirchhojJ energy 


117 


Note the similarity between this expression and the Saint- 




2 


MlldevaErf + ^tr^)] 2 , 


(2.7) 


where Ei = f (C — 1) = \{U 2 — 1) is the Green-Lagrangian strain. 

The squared Euclidean distance of F to SO(n) is often used as a lower bound for more general elastic 
energy potentials. Friesecke, James and Muller 78|, for example, show that if there exists a constant C > 0 
such that 

W(F) > C ■ distE Ucli d(r 1 , SO(3)) (2.8) 

for all F £ GL + (3) in a large neighbourhood of 1, then the elastic energy W shows some desirable properties 
which do not otherwise depend on the specific form of W. As a starting point for nonlinear theories of 
bending plates, Friesecke et al. also use the weighted squared norm 

II VWF- t\\l^ K = u. ||dev 3 (J7 - 1)|| 2 + f [tr(C7 - l)] 2 = H \\U - if + ^ [tr(t7 - l)] 2 , 


9 














where A is the first Lame parameter, as an energy function satisfying (|2.8|). The same energy, also called the 
Biot energy 


149 , has been recently motivated by applications in digital geometry processing [43 


However, the resulting strain measure u(U) = distEuciid(-F, SO(n)) = ||Z7 — 1|| does not truly seem appro¬ 
priate for finite elasticity theory: for U —> 0 we find \\U — 1|| —► ||1|| = ^/n < oo, thus singular deformations do 
not necessarily correspond to an infinite measure ui. Furthermore, the above computations are not compatible 

in general 149, 70, 711, 


with the weighted norm introduced in Section 2.1 


mm 

QeSO(ri) 


min || Q 
QeSO(n) 


F -111 ^ 


WF T F~ : 


I fl,g c ,K, ’ 


(2.9) 


thus the Euclidean distance of F to SO(ra) with respect to || . || K does not equal \\VF T F — 1|| K in 
general. In these cases, the element of best approximation is not the orthogonal polar factor R = polar(F). 


In fact, the expression on the left-hand side of (2.9) is not even well defined in terms of linear mappings 
F and Q 149 : the deformation gradient F = \7ip at a point x £ fl is a two-point tensor and hence, in 


particular, a linear mapping between the tangent spaces T x £l and T v i x \ip(Vl). Since taking the norm 


11*1 


= p ||dev n symX|| 2 + p c ||skew A|| 2 


[tr(X)] 2 


of X requires the decomposition of A' into its symmetric and its skew symmetric part, it is only well defined 
if X is an endomorphism on a single linear space J^] Therefore \\F ~ Q\\^ ^ K , while being a valid expression 
for arbitrary matrices F, Q £ R nxra , is not an admissible term in the setting of finite elasticity. 

We also observe that the Euclidean distance is not an intrinsic distance measure on GL + (n): in general, 
A- B $ GL+(n) for A,B £ GL+ (n), hence the term ||A — B\\ depends on the underlying linear structure 
of R nxn . Since it is not a closed subset of K nxn , GL + (n) is also not complete with respect to distEuciid; for 
example, the sequence (’ • l) N is a Cauchy sequence which does not converge. 

Most importantly, because GL + (n) is not convex, the straight line {A + t(B — A) \ t £ [0,1]} connecting 
A and B is not necessarily contained^] in GL + (n), which shows that the characterization of the Euclidean 
distance as the length of a shortest connecting curve is also not possible in a way intrinsic to GL + (n), as the 
intuitive sketcheffjin Figures [4] and [S] indicate. 

These issues amply demonstrate that the Euclidean distance can only be regarded as an extrinsic distance 
measure on the general linear group. We therefore need to expand our view to allow for a more appropriate, 
truly intrinsic distance measure on GL + (n). 


3 The Riemannian strain measure in nonlinear isotropic elasticity 

3.1 GL + (n) as a Riemannian manifold 

In order to find an intrinsic distance function on GL + (n) that alleviates the drawbacks of the Euclidean 
distance, we endow GL(n) with a Riemannian metric J^] Such a metric g is defined by an inner product 

g A : T a GL(rc) x T A GL(ra) -» R 

14 If X : V\ —> V 2 is a mapping between two different linear spaces Vi,V 2 , then X T is a mapping from V 2 to Vi, hence 
symX = | (X + X T ) is not well-defined. 

15 The straight line connecting F £ GL + (n) to its orthogonal polar factor R (i.e. the shortest connecting line from F to 
SO(n)), however, lies in GL"*" (n), which easily follows from the convexity of Sym +(n): for all t E [0,1], t U + (1 — t) 1 £ Sym + (n) 
and thus 

R + t{F - R) = R (t U + (1 - t) 1) £ R • Sym+(n) C GL+(n). 

16 Note that the representation of GL - *" (n) as a sphere only serves to visualize the curved nature of the manifold and that 
further geometric properties of GL”*" (n) should not be inferred from the figures. In particular, GL + (n) is not compact and the 
geodesics are generally not closed. 

17 For technical reasons, we define g on all of GL(n) instead of its connected component GL + (n); for more details, we refer to 
[129] , where a more thorough introduction to geodesics on GL(n) can be found. Of course, our strain measure depends only on 
the restriction of g to GL”*” (n). 
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Figure 4: The Euclidean distance as an extrinsic measure on GL' (n). 


on each tangent space T 4 GL (n), A £ GL(?t). Then the length of a sufficiently smooth curve 7 : [0,1] —> GL(rc) 
is given by 

L (i) = \J Wt) OH*), 7(f)) d t > 

where 7 (t) = ^ 7 (t), and the geodesic distance (cf. Figure I5J) between A,B £ GL + (n) is defined as the 
infimum over the lengths of all (twice continuously differentiable) curves connecting A to B: 

distgeod(A-B) = inf{L( 7 ) |7 £ C 2 ([0, l];GL + (n)), 7 (0) = A, 7 (1) = Bj . 

Our search for an appropriate strain measure is thereby reduced to the task of finding an appropriate 

GL+(n) / ad 4 


Figure 5: The geodesic (intrinsic) distance compared to the Euclidean (extrinsic) distance. 

Riemannian metric on GL(n). Although it might appear as an obvious choice, the metric g with 

g A {X, Y) := {X, Y) for all A £ GL+(n), X, Y £ R” xn (3.1) 

provides no improvement over the already discussed Euclidean distance on GL + (n): since the length of a 
curve 7 with respect to g is the classical (Euclidean) length 



- distg UC jj d (^4, B) = \\A — B\\ 2 

■- di Stgeod(X B ) 


L h} = J o \/9At) OW. i(t)) dt = J o IhWII dt ’ 

the shortest connecting curves with respect to g are straight lines of the form t i —> A + t{B — A) with 
A,B £ GL+ (n). Locally, the geodesic distance induced by g is therefore equal to the Euclidean distance. 
However, as discussed in the previous section, not all straight lines connecting arbitrary A,B£ GL + (n) are 
contained within GL + (n), thus length minimizing curves with respect to g do not necessarily exist (cf. Figure 
[ 6 ]). Many of the shortcomings of the Euclidean distance therefore apply to the geodesic distance induced by 
g as well. 

In order to find a more viable Riemannian metric g on GL(n), we consider the mechanical interpretation 
of the induced geodesic distance dist geo( j: while our focus lies on the strain measure induced by g , that is 
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Figure 6 : The shortest connecting ( geodesic ) curves in GL+(n) with respect to the Euclidean metric are straight lines, thus not 
every pair A, B £ G Jy (n) can be connected by curves of minimal length. The length of the straight line 7 : t 1—► A + 1 ,(B — A) 
connecting A to B is given by fg -\J 9-y{t) ( 7 (f), 7 (f)) dt = I B A | . whereas the curve 7 connecting A to C is not contained 
in GL+(n); its length is therefore not well defined. 


the geodesic distance of the deformation gradient F to the special orthogonal group SO (n), the distance 
dist g eod(-F| • F 2 ) between two deformation gradients F\,F 2 can also be motivated directly as a measure of 
difference between two linear (or homogeneous ) deformations F\, F 2 of the same body Q. More generally, we 
can define a difference measure between two inhomogeneous deformations ipi , <p 2 : C K n — > K n via 


dist(^i, <p 2 ) ■= / distgeod(V</Ji(x), Xip 2 {x))dx 

Jo. 


(3.2) 


under suitable regularity conditions for tpi,g >2 (e.g. if are sufficiently smooth with detVyjj > 0 up to 

the boundary). This extension of the distance to inhomogeneous deformations is visualized in Figure [7] 

In order to find an appropriate Riemannian metric g on GL(n), we must discuss the required properties 
of this “difference measure”. First, the requirements of objectivity (left-invariance) and isotropy (right- 
invariance) suggest that the metric g should be bi-0(n)-invarianf i.e. satisfy 

isotropy 

9 QA{QX, QY) = g A (X , Y) = 7aq(XQ , YQ) (3.3) 

' s * V* 

objectivity 

for all Q € O(n), A £ GL(n) and X, Y £ TaGL(u), to ensure that dist geo d(^4, B) = dist geo d (QA,QB) = 
distg eo d (^4. Q , B Cf ). 

However, these requirements do not sufficiently determine a specific Riemannian metric. For example, 
( |3.3| ) is satisfied by the metric g defined in (£3 as well as by the metric g with g 4 (X, Y) = ( A T X , A T Y). In 
order to rule out unsuitable metrics, we need to impose further restrictions on g. If we consider the distance 
measure dist(<^i, ipf) between two deformations introduced in (3.2), a number of further invariances 

can be motivated: if we require that the distance is not changed by the superposition of a homogeneous 
deformation, i.e. that 

dist(R • ipi,B ■ ip 2 ) = dist(</?i, v? 2 ) 

for all constant B £ GL(n), then g must be left-GLi(n)-invariant, i.e. 

g BA (BX,BY) = g A (X,Y) (3.4) 

for all A,B £ GL(n) and X, Y £ T A GL(n). The physical interpretation of this invariance requirement is 
readily visualized in Figure [ 8 ) 


It can easily be shown 


129 


invariant if and only if g is of the form 

g A (X,Y) = (A- 1 X,A~ 1 Y) 


that a Riemannian metric g is left-GL(n)-invariantpq as well as right-O(n)- 

«, (3-5) 


18 Of course, the left-GL(r).)-invariance of a metric also implies the left-0(n)-invariance. 
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Figure 7: The distance dist(y?i, (£ 2 ) : = /q dist geo d(V</?i(:r), Vcp 2 (x)) da: measures how much two deformations (pi,<p 2 of a body 
differ from each other via integration over the pointwise geodesic distances between Xipi(x) and 


where is the fixed inner product on the tangent space gl(ro) = T± GL(ra) = M" x " at the identity 

with 


( X , Y)^ tfJj0}K = 11 (dev n sym X, dev n sym Y) + g, c (skew X, skew Y) + | tr(X) tr(Y) 


(3.6) 


for constant positive parameters /x,/k c ,k > 0, and where (X,Y) = tr (X T Y) denotes the canonical inner 
product on gl(n) = K nx "0 A Riemannian metric g defined in this way behaves in the same way on all 
tangent spaces: for every A G GL + (n), g transforms the tangent space Ta GL + (n) at A to the tangent space 
T± GL + (n) = g((n) at the identity via the left-hand multiplication with A -1 and applies the fixed inner 
product (•, -)n^ 0 , K on fll(n) to the transformed tangents, cf. Figure [9] 

In the following, we will always assume that GL(n) is endowed with a Riemannian metric of the form 
(3.51 unless indicated otherwise. 

In order to find the geodesic distance 


dist g eod(-F, SO(n)) 


inf distgeod (A, 

QeSO(n) B 


of F G GL + (n) to SO(rc), we need to consider the geodesic curves on GL + (n). It has been shown 129 134 
5] that every geodesic on GL + (n) with respect to the left-GL(n)-invariant Riemannian metric induced by 


87 


the inner product (3.6) is of the form 


7 p(t) = F exp(t(sym£— ^skew£)) exp(f(l + ^)skew£) (3.7) 

19 If fi = fie = 1 and k = then the inner product (•, •)^,/x c ,k is the canonical inner product, and the corresponding metric 
g is the canonical left-invariant metric on GL(n) with gj\(X,Yf) = (A -1 X, A~ 1 Y) = tr (X T A~ T A~ 1 Y). Note that this metric 
differs from the trace metric gA(X,Y) = tr(A — 1 XA~ 1 T), cf. [551. 
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Figure 8 : The distance between two deformations should not be changed 
transformation B : dist(</?i, ^> 2 ) = dist (B • ipi, B • ^> 2 )- 


by the composition with an additional homogeneous 


with F £ GL + (n) and some £ £ gl(n), where exp denotes the matrix exponential p°| These curves are 
characterized by the geodesic equation 


c = 


P+Pc 

2p> 


(C r C-CC T ), 


C 7 l i- 


(3.8) 


Since the geodesic curves are defined globally, GL + (re) is geo desically complete with respect to the metric g. 
We can therefore apply the Hopf-Rinow theorem 111 129 to find that for all P, P £ GL + (n) there exists 
a length minimizing geodesic jp connecting F and P. Without loss of generality, we can assume that is 
defined on the interval [0,1]. Then the end points of jp are 


7f(0) = F and P = 7^(1) = F exp(sym£ — ^ skew£) exp((l + skew£), 


and the length of the geodesic 7 ^ starting in F with initial tangent F £ £ TpGL + (n) (cf. (3.71 and Figure 


11) is given by 129 


Lhi) = ll£ll. 


The geodesic distance between F and P can therefore be characterized as 

dist geod (F,P) = min{||£|| M)MciK | £ £ fl[(n) : 7>(1) = P} , 
that is the minimum of ||£|| M K over all £ £ fll(n) which connect F and P, i.e. satisfy 


exp(sym£ — skew£) exp((l + skew£) = F l P . 

lL 


(3.9) 


Although some numerical computations have been employed [21| to approximate the geodesic distance 
in the special case of the canonical left-GL(n)-invariant metric, i.e. for p = /j, c = 1, k = ^, there is no known 
closed form solution to the highly nonlinear system (3.9) in terms of £ for given F, P £ GL + (n) and thus no 
known method of directly computing dist geod (F, P) in the general case exists. However, this parametrization 
of the geodesic curves will still allow us to obtain a lower bound on the distance of F to SO(n). 


2t) The mapping £ 1 —>■ exp geod (£) := 7 ).( 1 ) = F exp(sym£ — 11,1 skew£) exp((l + IJc ) skew£) is also known as the geodesic 
exponential function at F. Note that in general exp geod (£) y F • exp(£) if £ is not normal (i.e. if 7 G C)- thus the geodesic 
curves are generally not one-parameter group s of t he form t 1 —>■ F exp(f £), hi contrast to bi-invariant metrics on Lie groups (e.g. 
SO(n) with the canonical bi-invariant metric 11361). 
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9a{X,Y) = (A-^A-'Y)^" 



Figure 9: A left-GL(n)-invariant Riemannian metric on GL(n) transforms the tangent space at A E GL”*~ (n) to the tangent 
space T± GIF (n) = 0 ((n) at the identity and applies a fixed inner product on g[(n) to the transformed tangents. Thus no 
tangent space is treated preferentially. 


3.2 The geodesic distance to SO(n) 

Having defined the geodesic distance on GL + (n), we can now consider the geodesic strain measure, which is 
the geodesic distance of the deformation gradient F to SO(ra): 


dist geo d(-F, SO(n)) — inf distg eoc j(T, Cf '). 

QeSO(n) 


(3.10) 


Without explicit computation of this distance, the left-GL(n)-invariance and the right-0(n)-invariance of 
the metric g immediately allow us to show the inverse deformation symmetry of the geodesic strain measure: 


dist g eod(T, SO(rc)) = inf dist geo d (F,Q) 

QeSO(n) 


= inf distgeod F,F 1 Q) 

QeSO(n) 


= inf dist g eod(l,T Q) = inf dist ge od(<3 Q,-F Q) 

QGSO(n) 8 QeSO(n) 8 


= inf distgeod(Q T ,-F X ) = distg eo d(-F \ SO (n)). 
QeSO(n) 


(3.11) 


This symmetry property demonstrates at once that the Eulerian (spatial) and the Lagrangian (referential) 
points of view are equivalent with respect to the geodesic strain measure: in the Eulerian setting, the inverse 
E -1 of the deformation gradient appears more naturalhp^] whereas F is used in the Lagrangian frame (cf. 
Figure 10). Equality (3.11) shows that both points of view can equivalently be taken if the geodesic strain 
measure is used. As we will see later on (Remark |3.5[ ), the equality dist geo d(H, SO(n)) = dist ge od(C', SO(n)) 
also holds for the right Cauchy-Green deformation tensor C = F T F = U 2 and the Finger tensor B = FF T = 
V 2 , further indicating the independence of the geodesic strain measure from the chosen frame of reference. 
This property is, however, not unique to geodesic (or logarithmic) strain measures; for example, the Frobenius 
norm 


\\E l/2 (U)\\ = \\\U-U- 




of the Bazant approximation Ey 2 = \ (U — U 1 ), cf. (1.4) and 17 , which can be considered a “quasiloga- 
rithmic” strain measure, fulfils the inverse deformation symmetry as wellrd However, it is not satisfied for 


Note that Cauchy originally introduced the tensors C 1 and B 1 in hi s invest igations of the nonlinear strain [4l[ |42[ fn\ 


2 The quantity -^||C7 — U 


|l82| , where C = F T F = U 2 is the right Cauchy-Green deformation ten sor 
Piola also formulated an early nonlinear elastic law in terms of C ~ 1 

1 11TT — rr 111 is suggested as a measure of strain magnitude by Truesdell and Toupin 204 p. 266]. 


8l][77] and B = FF T = V 2 is the Finger tensor, 
p. 347]. 
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Lagrangian frame 
(material setting) 




Eulerian frame 
(spatial setting) 


Figure 10: The Lagrangian and the Eulerian point of view are equivalently represented by the geodesic strain measure: 
distgeodl^l SO(n)) = dist geod (F -1 ,SO(n)). 


the Euclidean distance to SO(n): in general, 


\U — 1|| = distEuciid (-F, SO(n)) ^ dist Euclid (7 1 1 , SO(n)) = \\V 1 - 1|| . 


(3-12) 


Now, let F = RU denote the polar decomposition of F with U £ Sym + (n) and R £ SO(n). In order to 
establish a simple upper bound on the geodesic distance dist geo d(I 7 ', SO(n)), we construct a particular curve 
connecting F to its orthogonal factor R £ SO (n) and compute its length L( 7 #). For 

7 R(t) ■= R exp((l - t) log U ), 

where log U £ Sym(n) is the principal matrix logarithm of U, we find 

7fl(0) = R exp(log!7) = RU = F and 7 ^( 1 ) = R exp(0) = R £ SO(n). 


It is easy to confirm that 7 # is in fact a geodesic as given in (3.7) with £ = log U £ Sym(n). Since 
7fl 1 (*)7ii(*) = (R exp((l - t) log U))^ 1 R exp((l - t) log U) • (-log 17) = - log 17, 
the length of 7 # is given by 

l (1r) = J \J d* 

= \/(lR{t)- 1 iR(t),jR{t)~ lA fR(t))^c,n d t 

= J o \/(- 1 °S U, - log U)^ K dt = j 11 log 17|| M Atc K dt = 11 log f7|| MiMcjft . 

We can thereby establish the upper bound 


dist geod (F,SO(n)) = inf dist geod (F, Q) < dist geod (E, R) 

< L 2 (ir) = l|logt/||^ c)K = P ||dev„ log 17|| 5 


2 [tr(log 17)] 2 


(3.13) 


(3.14) 

(3.15) 


for the geodesic distance of F to SO(n). 

























Our task in the remainder of this section is to show that the right hand side of inequality (3.15) is also a 
lower bound for the (squared) geodesic strain measure, i.e. that, altogether, 


dist geod(^ S0 H) = H|dev„ log £/||" 


2 [tr(log U)f 


However, while the orthogonal polar factor R is the element of best approximation in the Euclidean case 
(for fi = fi c = 1, k = —) due to Grioli’s Theorem, it is not clear whether R is indeed the element in SO(n) 
with the shortest geodesic distance to F (and thus whether equality holds in (3.14)). Furthermore, it is not 
even immediately obvious that the geodesic distance between F and R is actually given by the right hand 
side of (3.15), since a shorter connecting geodesic might exist (and hence inequality might hold in ( 3.15| )). 

Nonetheless, the following fundamental logarithmic minimization property® of the orthogonal polar fac¬ 
tor, combined with the computations in Section 3.1 allows us to show that (3.15) is indeed also a lower bound 
for distgeod(-F)SO(n)). 


Proposition 3.1. Let F = R V F T F be the polar decomposition of F £ GL + (n) with R £ SO(n) and let 
denote the Frobenius norm on R nx ". Then 

inf ||symLog(Q T F)|| = ||symlog(i? T F)|| = ||log , 

QeSO(n) 


where 


n St Jl s y mL °g(<3 TF )ll : = _ inf inf{||symX| | X G R” xn , exp(X) =Q T F} 

QeSO(n) QeSO(n) 


is defined as the infimum of ||sym . || over “all real matrix logarithms’’ of Q T F. 


Proposition 3.1 which can be seen as the natural logarithmic analogue of Grioli’s Theorem (cf. Section 
2.2), was first shown for dimensions n = 2, 3 by Neff et al. 156 using the so-called sum-of-squared-logarithms 


inequality 29 170 48, 30 . A generalization to all unitarily invariant norms and complex logarithms for 
arbitrary dimension was given by Lankeit, Neff and Nakatsukasa 118 . We also require the following corollary 
involving the weighted Frobenius norm, which is not orthogonally invariant j 2 ^] 

Corollary 3.2. Let 


\\ x \\p,p c , K = HldeVnSymXH + p c ||skew X\Y 


MX )] 2 




for all X £ K" xn , where || . || is the Frobenius matrix norm. Then 

inf ||symLog(Q T F)|| = ||log VF T F\\ 

QGSO(n) 


p,g c ,K ' 


Proof. We first note that the equality detexp(A') = e tdAd holds for all X £ R” xn . Since det Q = 1 for all 
Q G SO(n), this implies that for all X £ M" xn with exp(X) = Q T F, 

tr(symX) = tr(A') = ln(det(exp(A))) = ln(det(Q T F)) = ln(detF). 

Therefor® 

||sym ^ ^ = n ||dev„ symX|| 2 + ^ [tr(symX)] 2 

= /r||symA|| 2 + n K ^ [tr(symA)] 2 = /i||symX|| 2 + U K ^ ^ (ln(detF)) 2 

23 Of course, the application of such minimization properties to elasticity theory has a long tradition: Leonhard Euler, in 
the appendix “De curvis elasticis” to his 1744 book “Methodus inveniendi lineas curvas maximi minimive proprietate gaudentes 
sive solutio problematis isoperimetrici latissimo sensu accepti” |62| |164| , already proclaimed that “[...] since the fabric of the 
universe is most perfect, and is the work of a most wise creator, nothing whatsoever takes place in the universe in which some 
rule of maximum and minimum does not appear 

24 While \\Q T XQ\\ K = ||X|| K for all X G E nXn and Q G O(n), the orthogonal invariance requires the equalities 
\\Q X \\= II X( 3II^, Mc , k = 1I- Y IU, Mc , k . which do not hold in general. 

25 Observe that p. ||dev n >'|| 2 + f [tr(F )] 2 = p ||F || 2 + n [tr(y )] 2 for all Y £ R nXn . 
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and finally 


QeSO(ra) \\ symLog ( QT ( 3 - 16 ) 
= inf inf{||symX||^ | X € M" x ” ) ex p(X) = q t f} 

QeSO(n) nrc, 

= inf inf{/i ||symX|| 2 + " K ~ 2/ " (ln(detF)) 2 I X G K nxn , exp(X) = Q T F} 

Qeso(n) 2 n 

= p inf inf{||symX|| 2 | X G K nxn , exp(X) = Q T F} + n ^ ~ (ln(detU)) 2 
QGSO(n) 2 Tl 

= M ||log VWf \\ 2 + (ln(detF)) 2 

= MlHog Vf^F \\ 2 + nK 2n 2fl [tr(logV^)] 2 

= T ||dev„ log VWF \\ 2 + f [tr(log V^)] 2 = ||log VWf\\^ k . ■ 

Note that Corollary |3.2| also implies the slightly weaker statement 

= Illog'^'lU,. 

by using the simple estimate ||X||^ c K > ||symX||^ c K . 


We are now ready to prove our main result. 

Theorem 3.3. Let g be the left-GL(n)-invariant, right-0(n)-invariant Riemannian metric on GL(n) defined 
by 

gA(x, y ) = (A _ 1 x, A~ 1 y) MiMe)K , h,h c ,k>q, 

for A G GL(n) and X , Y G K raxn , where 

(X, Y)^,* = p (dev n symX, dev n sym Y) + /r c (skewX, skew Y) + f tr(X) tr(Y). (3.17) 

Then for all F G GL + (n), the geodesic distance of F to the special orthogonal group SO (n) induced by g is 
given by 

distg eod (F, SO(n)) = p ||dev n log U\\ 2 + | [tr(log U )] 2 , (3.18) 

where log is the principal matrix logarithm, tr(X) = X,; d denotes the trace and dev„ X = X — I tr(X) • 1 

is the n-dimensional deviatoric part of X G R" x ”. The orthogonal factor R G SO(n) of the polar decomposi¬ 
tion F = RU is the unique element of best approximation in SO(n), i.e. 

distgeod(-F) SO(n)) = dist geo d(-F, R) = dist geo d (R T F, 1 ) = dist geo d(£7,1) • 


In particular, the geodesic distance does not depend on the spin modulus p c . 

Remark 3.4 (Uniqueness of the metric). We remark once more that the Riemannian metric considered in 
Theorem |3.3| is not chosen arbitrarily: every left-GL(n)-invariant, right-0(?z)-invariant Riemannian metric 
on GL(n) is of the form given in (3.17) for some choice of parameters p,p c ,n > 0 129 . □ 


Remark 3.5. Since the weighted Frobenius norm on the right hand side of equation (3.18) only depends 
on the eigenvalues of U = y/F T F, the result can also be expressed in terms of the left Biot-stretch tensor 
V = VFF T , which has the same eigenvalues as U: 


distgeod^ SO(n)) = p ||dev n log W|| 2 


-[tr(logU)] 2 


(3.19) 
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Applying the above formula to the case F = P with P £ Sym + (n), we find VP T P = \/PP T = P and 
therefore 

dist 2 (P, SO(n)) = dist 2 (P, 1) = /x ||dev„ log P|| 2 + ^ [tr(log P)] 2 , (3.20) 

since 1 is the orthogonal polar factor of P. For the tensors U and V, the right Cauchy-Green deformation 
tensor C = F T F = U 2 and the Finger tensor B = FF T = V 2 , we thereby obtain the equalities 

dist geo d(P, SO(rc)) = dist g eod(P, 1) = distg eo d(P _1 ,1) (3-21) 

— distgeod(C, 1) = distgeod(0 ,1) — distg eo d(G, SO(ra)) 
and dist g eod(V', SO(n)) = distgeod) 1 ^, 1) = dist geo d(^ _1 ,l) (3.22) 

= dist g eod(P? H) = dist g eod(b^ , 1) = distg eo d(P; SO(n)) . 



exp(sym£ — ^ skew£) exp((l + ^) skew£) = F 1 Q (3.23) 


and 

= distgeod(P) Q) ■ (3.24) 

In order to find a lower estimate on ||£|| K (and thus on dist ge od(P, Q)), we compute 


exp(sym£ — ^ skew£) exp((l + ^) skew£) = F 


l Q 


exp((l + ^) skew£) 


-l 


exp(sym£ - ^ skew£) 1 = Q 1 F 

exp(— sym£ + ^ skew£) = exp( (1 + ^) skew£) Q T F . 

"-V-" 

eso(n) 


Since exp(W) £ SO (n) for all skew symmetric W £ so(n), we find 



Figure 11: The geodesic (intrinsic) distance to SO(n); neither the element 
F £ G Tp GL+ (n) of the connecting geodesic is known beforehand. 


Q of best approximation nor the initial tangent 
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(3.25) 


exp(— sym£ + ^skew£) = QjF 


=:Y 


with = Q exp (—(1 + ^)skew£) € SO(?z); note that symF = — sym£. According to (3.25), Y = 
— sym£ -f — skew£ is “a logarithm’! 5 ®] of QjF. The weighted Frobenius norm of the symmetric part of 
F = - sym£ + ^ skew£ is therefore bounded below by the infimum of HsymXH^ ^ K over “all logarithms” 
X of QjF: 


IlsymCII 


,K 


= ||symF|| 


/X,/I C ,K 


Qe so( 

We can now apply Corollary |3.2| to find 

,2 /rr. II All 2 


inf{||symA|| MAtcif 

> inf inf{||sym X|| 

“ QeSO(n) 111 J 

= inf JsymLog(Q T F)|| 


\X G . 
\X G 


', exp(X) = QjF} 
' • exp(X) = Q t F} 


(3.26) 


dist geod (F,Q) = ll£ll MlMc! « = /z||dev n sym£|| + /z c ||skew£|| + - [tr(sym £)] 2 

> /z ||dev„ sym £|| 2 + ^ [tr(sym £)] 2 


(3.27) 


= ||sym£|| 




J3.261 


A inf ||symLog(Q J F)|| 


QeSO(n) 


Corollary |3.2| 


M 11 log ^F t F\\ 1^ cK = fi ||dev„ log 17|| 2 + * [tr(log U )] 2 


for U = \JF T F. Since this inequality is independent of Q and holds for all Q G SO (n), we obtain the desired 
lower bound 

distg eod (F, SO(n)) = ^ inf dist 2 eod (F, Q) > /z l|dev„ log U || 2 + ^ [tr(log U)] 2 
QeSO(n) 2 

on the geodesic distance of F to SO (n). Together with the upper bound 

distg eod (F,SO(n)) < distg eod (F, R) < \i ||dev„ log U\\ 2 + ^ [tr(logC /)] 2 


already established in (3.15), we finally find 

dist 2 eod (y, SO(n)) = dist 2 eod (F, R) = /z ||dev n log U\\ 2 


[tr(log U)] 2 . 


(3.28) 


By equation (3.28), apart from computing the geodesic distance of F to SO(n), we have shown that 
the orthogonal polar factor R = polar(F) is an element of best approximation to F in SO(?z). However, 
it is not yet clear whether there exists another element of best approximation, i.e. whether there is a Q G 
SO (n) with Q ^ R and dist geod (F, Q) = dist geod (F, R) = dist geo d(A, SO(n)). For this purpose, we need 
to compare geodesic distances corresponding to different parameters /z, /z c , k. We therefore introduce the 
following notation: for fixed /z, /z c , k > 0, let dist geo d,/r,Ai c ,K denote the geodesic distance on GL + (n) induced by 
the left-GL(n)-invariant, right-0 (?z)-invariant Riemannian metric g (as introduced in (3.5)) with parameters 
/z, fjL c , k. Furthermore, the length of a curve 7 with respect to this metric will be denoted by £ MiAtcjK ( 7 ). 


26 Loosely speaking, we use the term “a logarithm of A S GLd(n)” to denote any (real) solution X of the matrix equation 
exp X = A. 
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Assume that Q £ SO(n) is an element of best approximation to F with respect to g for some fixed 
parameters /i,p, c ,ze > 0. Then there exists a length minimizing geodesic 7 : [0,1] —► GL + (n) connecting Q to 
F of the form 

7 (t) = Q exp(f(sym£ — skew£)) exp(f(l + skew£) 


and the length of 7 is given by 

r 2 /.A ll^l|2 


with f e R nx ", 

L 1,U. C , K ( 7 ) = ll^l^ Mo , K = M ||dev n sym£|| 2 + /z c ||skew£|| 2 + ^ [tr(£)] 2 
We first assume that skew£ ^ 0. We choose /z c > 0 with jx c < fj, c and find 

distg eod ^ AciK (E,SO(n)) = inf dist geod,/ i ,/I c ,«( F ’ < 3) 

< dist 2 eod ^ jK (F,Q) < T^ Ao , k (7), 


(3.29) 


since 7 is a curve connecting F to Q £ SO(n); note that although 7 is a shortest connecting geodesic with 
respect to parameters /z, /z c ,k by assumption, it must not necessarily be a length minimizing curve with 

respect to parameters /x, fl c , k. Obviously, ||£|| - c K < ||£|| p K if skew£ ^ 0, and therefore 

LU 0 ,M = lien l,a c , K < lieilLc,« = = dist 2 eodj ^ cjK (F, Q ). (3.30) 

By assumption, Q is an element of best approximation to F in SO(n) for parameters /i, /z c , k, thus 

dist leod,^ c , K ( F ^ Q) = dist geod,/x,Mc,tc( F ; SO(n)) (3.31) 

= /z ||dev„ log U \\ 2 + | [tr(log U )} 2 = dist 2 eodjM>AcjK (F, SO(ra)), 

where the last equality utilizes the fact that the distance from F to SO (n) is independent of the second 

parameter (// c or fi c ). Combining (3.29), (3.30) and (3.31), we thereby obtain the contradiction 

distg eodiMiAcjK (F,SO(n)) < L 2 iAcjK ( 7 ) < dist 2 eodj/ij/icjK (.F, Q) = dist 2 eodj/JjAcjK (.F, SO(n)), 

hence we must have skew£ = 0. But then 

7 ( 1 ) = Q exp(sym£— ^skew£) exp((l + ^)skewf) = Q exp(sym£), 

and since exp(sym^) £ Sym + (n), the uniqueness of the polar decomposition F = RU yields exp(sym£) = U 
and, finally, Q = R. ■ 

The fact that the orthogonal polar factor R = polar(F) is the unique element of best approximation to 
F in SO(rz) with respect to the geodesic distance corresponds directly to the linear case (cf. equality (2.3) in 
Section |2.1[ ), where the skew symmetric part skew Vm of the displacement gradient Vzt is the element of best 
approximation with respect to the Euclidean distance: for F = 1 + Vzt we have 

U = t + sym Vzz + 0(||Vzz|| 2 ) and R = t + skew Vzz + 0(\\ Vzz|| 2 ), 

hence the linear approximation of the orthogonal and the positive definite factor in the polar decomposition 
is given by skew Vzz and sym Vzz, respectively. The geometric connection between the geodesic dis tan ce on 
GL + (zz) and the Euclidean distance on the tangent space R nxn = g[(n) at 1 is illustrated in Figure 


12 


Remark 3.6. Using a similar proof, exactly the same result can be shown for the geodesic distance dist geo d,right 


induced by the right-GL(n)-invariant, left-0(n)-invariant Riemannian metric 207 

g T it\X,Y) = (XA-\YA~ 1 )^ 
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Figure 12: The isotropic Hencky energy of F measures the geodesic distance between F and SO(n). The linear Euclidean strain 
measure is obtained as the linearization via the tangent space gl (n) at 1. 


on GL(n): 

dist g eod ,right (F, SO (n)) = distg eod (F, SO(ra)) = p ||dev„ log U\\ 2 + | [tr(log U )] 2 . 


The right-GL(n)-invariant Riemannian metric can be motivated in a way similar to the left-GL(n)-invariant 
case: it corresponds to the requirement that the distance between two deformations F\ and F 2 should not 
depend on the initial shape of fi, i.e. should not be changed if O is homogeneously deformed beforehand (cf. 
Figure [l3|. A similar independence from prior deformations (and so-called “pre-stresses”), called “elastic 
determinacy” by L. Prandtl 171 , was postulated by H. Hencky in the deduction of his elasticity model; cf. 

□ 


100, p. 618], 146, p. 19] and Section 4.2 


According to Theorem |3.3| the squared geodesic distance between F and SO (n) with respect to any 
left-GL(n)-invariant, right-0(n)-invariant Riemannian metric on GL(n) is the isotropic quadratic Hencky 
energy 


Wh(F) = p ||dev„ log U\[‘ 


[tr(log U)] 2 


where the parameters p, k > 0 represent the shear modulus and the bulk modulus, respectively. The Hencky 
energy function was introduced in 1929 by H. Hencky 101 , who derived it from geometrical considerations 


as well: his deductioiF] was based on a set of axioms including a law of superposition (cf. Section 4.2) for the 
stress response function 146 , an approach previously employed by G.F. Becker 18} 152] in 1893 and later 


followed in a more general context by H. Richter 176 , cf. 177 175 178 . A different constitutive model for 


uniaxial deformations based on logarithmic strain had previously been proposed by Imbert 114 and Hartig 


[89 . While Ludwik is often credited with the introduction of the uniaxial logarithmic strain, his ubiquitously 
cited article 124] (which is even referenced by Hencky himself [l02| p. 175]) does not provide a systematic 
introduction of such a strain measure. 


While the energy function Wn(F) = distg eod (A, SO(n)) already defines a measure of strain as described 
in Section [T t[ we are also interested in characterizing the two terms ||dev„ log U\\ and |tr(logf/)| as separate 
partial strain measures. 

27 Hencky’s approach is often misrepresented as empirically motivated. Truesdell claims that “Hencky himself does not give 
a systematic treatement ” in introducing the logarithmic strain tensor [199| p. 144] and attributes the axiomatic approach to 
Richter [176] instead j 204, p. 270]. Richter’s resulting deviatoric strain tensors dev 3 log U and rleva log V' are disqualified as 
“complicated algebraic functions” by Truesdell and Toupin [204] p. 270]. 
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Figure 13: The right-GL(n)-invariance of a distance measure on GL(n): the distance between two homogeneous deformations 
F\, F 2 is not changed by a prior homogeneous deformation B , i.e. dist geo d(Fi, F 2 ) = dist geo( j(-Fi • B, F 2 • B). 


Theorem 3.7 (Partial strain measures). Let 

w iso (F) := ||dev n log VF T F\\ and u vo i(F) := |tr(log VF T F) \ . 

Then 


w iso (F) = dist geodi SL(n) ( det F y n > SO(ra) 


and 


Wvoi(-F) = Vn ■ distgeod^+.i ((detF) 1/n ■ 1, l) , 

where the geodesic distances dist geodj sL(n) an d dist geodj r+.j on the Lie groups SL(n) = {A £ GL(n) | det A = 
1} and R + • 1 are induced by the canonical left-invariant metric 

9a(X,Y ) = {A~ 1 X 1 A~ l Y) = tr{X T A~ T A~ l Y). 


Remark 3.8. Theorem |3.7| states that Wj so and w vo i appear as natural measures of the isochoric and 
volumetric strain, respectively: if F = F; so F vo \ is decomposed multiplicatively 73 into an isochoric part 
Fiso = (det-F) -1 /" • F and a volumetric part F vo \ = (detF) 1 /™ • t, then Wj so (F) measures the SL(?z)-geodesic 
distance of Fi so to SO(n), whereas ~^uj vo \(F) gives the geodesic distance of F vo i to the identity 1 in the 
group R + • 1 of purely volumetric deformations. □ 


Proof. First, observe that the canonical left-invariant metrics on SL(n) and R + • 1 are obtained by choosing 
/i = /Lt c = 1 and k = ^ and restricting the corresponding metric g on GL + (n) to the submanifolds SL(n), 
R + • 1 and their respective tangent spaces. Then for this choice of parameters, every curve in SL(n) or R + • 1 
is a curve of equal length in GL + (n) with respect to g. Since the geodesic distance is defined as the infimal 
length of connecting curves, this immediately implies 


dist geodj sL(n) (Fiso, SO(n)) > dist geodjGL+(n) (F iso , SO(ra)) 
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as well as 


distg eoc j ) jj+.i (Fvoi, 1) ^ distg eoc j iG L+(n) (-Fvoh H) — disl'geod, ql+ ( F vo \ , SO (n)) 

to obtain the lower 


for F iso := (det F)~^ n ■ F and F v oi := (detF) 1 /" • 1. We can therefore use Theorem 
bound® 

dlstgeod, SL(n) (-^isoj SO(?7-)) 

> dist geod,GL+( n ) ( F iso, SO(fl)) 

= ||dev„ log (sjFT o F ls ^\ || 2 + i tr ("log y/F? 0 F is 


3.3 


and 


-i/ 


=i 


Hllog (det/^ 0 F iso ) y^F~ 0 || 2 + i In ( det yf^F, 


= ||log (y/l&f 1^) || 2 = 11log ((det F)~^VWf) || 2 = ujf so (F) 

dis tg eod ,R+.i (Fvoh 1) > dist 2 eod GL+( „ ) (F vo i, SO(rc)) 

= ||dev„ log ^Fj ol F voi y + * [tr(log (y/F^F m ^j )] 2 

= 11dev„ (ln((det F) 1/71 ) • l) || 2 + i [ln(det ((detF) 1 /" ■ l))] 5 


= - [In (det V F T F)] 2 = - [tr(log VF T F)] 2 = - w 2 , (F ). 
n n n 

To obtain an upper bound on the geodesic distances, we define the two curves 

7 iso: [0,1] SL(n), 7 iso (t) = R exp (t clev n log U) 


(3.32) 


(3.33) 


and 

7voi: [0,1] -► R+ ■ 1, 7 voi(t) = e " tr(1 ° sU) ■ 1, 

where F = RU with R £ SO(n) and U £ Sym + (n) is the polar decomposition of F. Then 7 ; so connects 
(det F) -1 /” • F to SO(n): 

7 iso( 0 ) = UeSO(n), 

7iso(l) = R exp(dev„ log U) = R exp(log U - tr(1 ° g u) ■ 1) 

= R exp (log U) exp(- tr(I ° sC/) • 1) 

= RU exp(-is^T . 1 ) = (det U)~^ n ■ F = (det F)~ 1/n ■ F, 
while 7 vo i connects (detF) 1 /" • 1 and 1: 

7vo i(0) = 1, 7voi(l)= e s tr(losC/) • 1 = e (r ln ( det u ) . 1 = (det U)^ n ■ 1 = (detF) 1 /" • 1. 


28 For some of the rules of computation employed here involving the matrix logarithm, we refer to LemmalA.llin the appendix. 
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The lengths of the curves compute to 


L{ 7 iso) = f || 7 iso(i) * 7 iso(i)|| dt ( 3 . 34 ) 

JO 

1 

|| (R exp {t dev„ log t/)) _1 R exp(f dev„ log U) dev n log ?7|| d t 


||dev„ log U || d t = ||dev n log VF T F\\ = w iso (F) 


as well as 

L( 7voi) = [ ||7voi(i) _1 7voi(i)|| dt (3.35) 

Jo 

||( e S tr (logl7) . t )-l . tr (logt/) . e n tr(log U) . t || dt 

= f • 111 d t = |tr(1 ° g[/)l • ||1|| = -L |tr(log Vf^F)\ = -L uj vo i(F) , 

J o n \ n v n 

showing that 




distgeod,SL(n) ((detF)-V" • F, SO (n)) < L 2 ( 7 is o) = < 4 G (F) 


and 


dist g eod ,H+-i ((detF) 1/n • 1, l) < L 2 ( 7 voi) = ^ • w 2 oi(F), 


which completes the proof. 

Remark 3.9. In addition to the isochoric (distortional) part -Fisc = (detF)^ 1 / 11 • F and the volumetric part 
F vol = (det F) 1 /” • 1, we may also consider the cofactor Cof F = (det F) ■ F T of F £ GL + (n). Theorem 
allows us to directly compute (cf. Appendix | A. 4|) the distance 


3.3 


distggod(Oof F, SO(n)) = g ||dev„ log C/|| 2 


k {ji — l ) 2 


[tr(log U)f 


□ 


4 Alternative motivations for the logarithmic strain 

4.1 Riemannian geometry applied to Sym+(n) 

Extensive work on the use of Lie group theory and differential geometry in continuum mechanics has already 
been done by Rou gee 181 180 |182 183|, Moakher 


[65j [66] |67| [68] (cf. |11911120] 116311167 


137 139 , Bhatia 26 


64 


symmetric (3 x 3)-tensors with the Riemannian metricr® 


and, more recently, by Fiala 
166 ). They all endowed the convex cone Sym + (3) of positive definite 


gc(X,Y) = ti(C- 1 XC~ 1 Y) = (XC~\C~ l Y) = {C~^ 2 X C~ 1/2 , C -1/2 Y C~ 1/2 ) 


(4.1) 


where C € Sym+(3) and A, Y € Sym(3) = TcSym+(3). Fiala and Rougee deduced a motivation of the 
logarithmic strain tensor log U via geodesic curves connecting elements of Sym + (?r). However, their approach 

29 Note the subtle difference with our metric gc(X. V) = (C~ 1 X,C~ 1 Y). Pennec |166| p. 368] generalizes i |4. 1 [ l by using the 
weighted inner product (X,Y), = (X,Y) + (3 tr(X) tr(Y) with f) > — L. 
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differs markedly from our method employed in the previous sections: the manifold Sym + (n) already corre¬ 
sponds to metric states C = F T F , whereas we consider the full set GL + (n) of deformation gradients F (cf. 
Appendix A.3 and Table [l] in Section [6]). This restriction can be viewed as the nonlinear analogue of the a 
priori restriction to e = syrnVit in the linear case, i.e. the nature of the strain measure is not deduced but 
postulated. Note also that the metric g cannot be obtained by restricting our left-GL(3)-invariant, right- 
0(3)-invariant metric g to Sym + (3)p°| Furthermore, while Fiala and Rougee aim to motivate the Hencky 
strain tensor log U directly, our focus lies on the strain measures Wi so , w vo i and the isotropic Hencky strain 
energy W H . 

The geodesic curves on Sym + (n) with respect to g are of the simple forn:p~| 


7(i) = cf exp(t • C~ 1 / 2 MC~ 1/2 ) cl /2 


(4.2) 


with C\ £ Sym + (n) and M £ Sym(n) = Tc 1 Sym + (n). These geodesics are defined globally, i.e. Sym + (n) 
is geodesically complete. Furthermore, for given C ±, C*2 £ Sym + (n ), there exists a unique geodesic curve 
connecting them; this easily follows from the representation formula ( |4.2[ ) or from the fact that the curvature 
of Sym + (n) with g is constant and negative 65 116 25 . Note that this implies that, in contrast to GL + (n) 


with our metric g, there are no closed geodesics on Sym + (n). 

An explicit formula for the corresponding geodesic distance was given by Moakherru 


dist geodl Sym+(n)(C'l ) C 2 ) = ||log(C 2 1/2 C X C. 




(4.3) 


In the special case C 2 = 1, this distan ce m easure is equal to our geodesic distance on GL + (n) induced by 
the canonical inner product: Theorem 3.3 applied with parameters g = fi c = 1 and k = j- to R = 1 and 
U = Ci, shows that 

distgeod, GL+ (n) (£l) -^-) I log \ \ distg eod) Sy m + (n) (^1, H) ■ 

More generally, assume that the two metric states C \, C -2 £ Sym + (n) commute. Then C ^ 1 C'i £ Sym + (n), 
and the left-GL(n)-invariance of the geodesic distance implies 


dist geodjG L+(n)(Cij C 2 ) — dist geod!GL +( n )(C 2 1 C'i,l) — ||log(C 2 1 C'i)|| 

= ||log(C 2 - 1/2 C 2 - 1/2 C 1 )|| = ||log (C~ 1/2 CiC- 1 ' 2 ) 
dist geo( p Sym+(n)(f^l; G2) . 


(4.4) 


However, since C^Ci ^ Sym + (n) in general, this equality does not hold on all of Sym + (n). 

A different approach towards distance functions on the set Sym + (n) was suggested by Arsigny et al. [8j[9, 
71 who, motivated by applications of geodesic and logarithmic distances in diffusion tensor imaging, directly 
define their Log-Euclidean metric on Sym + (n) by 


dist L og-Euciid(C'i,C' 2 ) := ||log C\ — logC^H 


(4.5) 


197 

with 


30 Since Sym + (n) is not a Lie group with respect to matrix multiplication, the metric g itself cannot be left- or right-invariant 
in any suitable sense. 

31 While Moakher gives the parametrization stated here, Rougee writes the geodesics in the form 7 (t) = exp(t-Log(C 2 C^~ )) C\ 

with Ci,C 2 £ Sym”*"(n), which can also be written as 7 (t) = )*Ci; a similar formulation is given by Tarantola 

eq. (2.78)]. For a suitable definition of a matrix logarithm Log on GL"*~(n), these representations are equivalent to ( |4.2| ) 

M = log (C ~ 1/2 Ci C“ 1/2 ) £ Sym(n). 

32 Moakher 137 eq. (2.9)] writes this result as ||Log(C^~ 1 Ci)|| = 1 1h 2 where A i are the eigenvalues of C^” 1 Ci. The 

right hand side of this equation is identical to the result stated in (|4.3|). However, since C7 1 Ci is not necessarily normal, there 


is in general no logarithm Log(C 2 L C\) whose Frobenius norm satisfies this equality. Note that the eigenvalues of the matrix 
C^Ci are real and positive due to its similarity to C 2 ^ 2 (C ^ -1 Ci)C 2 ^ 2 = C 2 1 ^ 2 CiC 2 ^ 2 £ Sym + (n). 
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where || . || is the Frobenius matrix norm. If G\ and C 2 commute, this distance equals the geodesic distance 
on GL + (n) as well: 

dist geod , GL+(n) (C' 1) C' 2 ) = ||log(C 2 - 1 C 1 )|| 

= ||log(C 2 - 1 ) +log^OH (4.6) 

= ||log Ci - log C 2 || = dist L og-Euciid(Ci,C 2 ), 


where equality in (4.6 1 holds due to the fact that C\ and C 2 commute. Again, this equality does not hold 
for arbitrary C\ and C 2 . 

Using a similar Riemannian metric, geodesic distance measures can also be applied to the set of positive 
definite symmetric fourth-order elasticity tensors, which can be identified with Sym + (6). Norris and Moakher 
applied such a distance function in order to find an isotropic elasticity tensor C: Sym(3) —> Sym(3) which 


best approximates a given anisotropic tensor 138 157 


The connection between geodesic distances on the metric states in Sym + (?r) and logarithmic distance 
measures was also investigated ex tensi vely by the late Albert Tarantola 197 , a lifelong advocate of logarithmic 
measures in physics. In his view 197 p. 4.3.1], “... the configuration space is the Lie group GL + (3), and the 


only possible measure of strain (as the geodesics of the space) is logarithmic.'” 


4.2 Further mechanical motivations for the quadratic isotropic Hencky model 
based on logarithmic strain tensors 

“At the foundation of all elastic theories lies the definition of strain, and before introducing a new 
law of elasticity we must explain how finite strain is to be measured .” 

Heinrich Hencky: The elastic behavior of vulcanized rubber [103]. 


Apart from the geometric considerations laid out in the previous sections, the Hencky strain tensor E 0 = log U 
can be characterized via a number of unique properties. 

For example, the Hencky strain is the only strain tensor (for a suitably narrow definition, cf. 152] ) that 
satisfies the law of superposition for coaxial deformations: 

Eo{Ui ■ U 2 ) = Eoilh) + E 0 (U 2 ) (4.7) 

for all coaxial stretches U\ and U 2 , i.e. U\,U 2 £ Sym + (n) such that U\ ■ U 2 = U 2 ■ U\. This characterization 


was used by Heinrich Hencky 196 97 102 103 in his original introduction of the logarithmic strain tensor 
99] 101, 100][l46j and, indeed much earlier, by the geologist George Ferdinand Becker 133 , who postulated 


a similar law of superposition in order to deduce a logarithmic constitutive law of nonlinear elasticity 18 
|152] (cf. Appendix |A.2[ ). 

In the case n = 1, this superposition principle simply amounts to the fact that the logarithm function 


/ = log satisfies Cauchy’s 40 well-known functional equation 


/(Ai • A 2 ) = f(Xi) + /(A 2 ), 


(4.8) 


i.e. that the logarithm is an isomorphism between the multiplicative group (R+, •) and the additive group 
(M, +)• This means that for a sequence of incremental one-dimensional deformations, the logarithmic strains 
ef og can be added in order to obtain the total logarithmic strain e[°* of the composed deformation 
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e iog 


■ e log 


• ■ • + efog = log y- + log y- + . . . + log 

-^0 Lr\. L'n-l 


= log ~r~ = efog, 

-L/f) 


where Li denotes the length of the (one-dimensional) body after the i-th elongation. This property uniquely 
characterizes the logarithmic strain ei og among all differentiable one-dimensional strain mappings e: R + —> R 
with e'(l) = 1. 
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Since purely volumetric deformations of the form A-l with A > 0 are coaxial to every stretch U £ Sym + (ro), 
the decomposition property (4.71 allows for a simple additive volumetric-isochoric split of the Hencky strain 


tensor 176 


log U = log 


U 


(det U) 1 /™ 


isochoric 


(det U) 1/n ■ 1 


= log 


U 


(det U)V n 


+ log 


(det U) 1/n ■ 1 


volumetric 


= dev n log V H-tr(log U) • t . 

n 


isochoric 


volumetric 


In particular, the incompressibility condition det F = 1 can be easily expressed as tr(log U) = 0 in terms of 
the logarithmic strain tensor. 

4.2.1 From Truesdell’s hypoelasticity to Hencky’s hyperelastic model 

As indicated in Section o the quadratic Hencky energy is also of great importance to the concept of 


hypoelasticity 


83 


Chapter IX]. It was found that the Truesdcll equatior 33 
[r] = 2 /x D + A tr(D) • 1 


199 


1201 


200 
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D = sym(F F x ), 


(4.9) 


with constant Lame coefficients /x, A > 0, under the assumption that the stress rate ^ 

rfl • , , .. • ... , . d l0 S 


IS 




and 


corotational , is satisfied if and only if 4_is the so-called logarithmic corotational rate ^ b and r = 2 /x log V+ 

A tr(logV)-l 211 209 158 172 173,212 213 214] , i.e. if and only if the hypoelastic model is exactly Hencky’s 
hyperelastic constitutive model. Here, r = det F-a{V) denotes the Kirchhoff stress tensor and D is the unique 
rate of stretching tensor (i.e. the symmetric part of the velocity gradient in the spatial setting). A rate f t 
is called corotational if it is of the special form 


d t 


[X] = X — fiX + XCl withS2eso(3) 


(4.10) 


which means that the rate is computed with respect to a frame that is rotated! 33 ] This extra rate of rotation 


is defined only by the underlying spins of the problem. Upon specialisation, for p = 1, A = 0 we obtai 1 ^ 33 
eq. 71] 

d log 


34 


s i l0 ^] = D 


as the unique solution to (4.9) with a corotational rate. 


Note that this characterization of the spatial 

p. 49, 


logarithmic strain tensor log V is by no means exceptional. For example, it is well known that 
Theorem 1.8] (cf. [35] ) 

d A 

— [A]=A + L T A + AL = D , 


90 


33 It is telling to see that equation ( |4.9| ) had already been proposed by Hencky himself in 100| for the Zaremba-Jaumann 
stress rate (cf. < |4.13| )). Hencky’s work, however, contains a typographical error [l00| eq. (10) and eq. (He)] changing the order 


of indices in his equations (cf. [33| ). The strong point of writing ( |4.9| ) is that no discussion of any suitable strain tensor is 
necessary. 

34 A rate is called objective if [S(QBQ T )] = Q (-^[S(B)])Q T for all (not necessarily constant) Q = Q(t) E O(n), 
where S is any objective stress tensor, and if =0^-5 = 0, i.e. th e motion is rigid if and only if = 0. 


35 Corotational rates are also special cases of Lie derivatives 


112 


127 


36 Cf. Xiao, Bruhns and Meyers [210 p. 90]: “.. . the logarithmic strain [does] possess certain intrinsic far-reaching properties 
[which] establish its favoured position in all possible strain measures ”. 
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where A = = |(1 — B 1 ) is the spatial Almansi strain tensor and 4t' is the upper Oldroyd rate (as 


defined in (4.14)). 


The quadratic Hencky model 

t = 2 fi log V + \ tr(log V) • 1 = AogvWnQog V) 
was generalized in Hill’s generalized linear elasticity laws ^ 108 eq. (2.69)] 


T r = 2 /x E r + A tr (E r ) • 1 


(4.11) 


(4.12) 


with work-conjugate pairs (T r ,E r ) based on the Lagrangian strain measures given in (1.3); cf. Appendix A.2 


for examples. The concept of work-conjugacy was introduced by Hill 106 
spatial stress power must be equal to its Lagrangian counterpart: 

det F ■ ( a , D) = (T r , E r ), 


via an invariance requirement; the 


(work-conjugacy) 


by means of which a material stress tensor is uniquely linked to its (material rate) conjugate strain tensor. 
Hence it generalizes the virtual work principle and is the foundation of derived methods like the finite element 
method. 

For the case of isotropic materials, Hill 106 p. 242] (cf. 109 ) shows by spectral decomposition techniques 
that the work-conjugate stress to log U is the back-rotated Cauchy stres s a multiplied by det F, hence (a, D) = 

additionally found 


106 


108 


185 


Sansour 

£>2 denotes the second Piola- 


tiere 


(. R 1 a R, ^ log U), which is a generalization of Hill’s earlier work 
that the Eshelby-like stress tensor S = CS 2 is equally conjugate toTogT7i" 

Kirchhoff stress tensor. For anisotropy, however, the conjugate stress exists but follows a much more complex 
format than for isotropy 109 . The logarithm of the left stretch log V in contrast exhibits a work conjugate 
stress tensor only for isotropic materials, namely the Kirchhoff stress tensor r = detF • a 162 l09 . 

While hyperelasticity in its potential format avoids rate equations, the use of stress rates (i.e. stress 
increments in time) may be useful for the description of inelastic material behavior at finite strains. Since 
the material time derivative of an Eulerian stress tensor is not objective, rates for a tensor X were developed, 
like the (objective and corotational) Zaremba-Jaumann rate 


d t 


[X] = X - WX + XW , W = skew L , L = FF~ 


(4.13) 


or the (objective but not corotational) lower and upper Oldroyd rates 


— [X] = X + L t X + XL 


and 


— \X] = X-LX- XL t . 
dt 1 1 


(4.14) 


to name but a few (cf. 90, Section 1.7] and 186]). Which one of these or the great number of other objective 
rates should be used seems to be rather a matter of taste, hence of arbitrariness] or heuristicSI but not a 
matter of theory. 

The concept of dual variables^ |as introduced by Tsakmakis and Haupt in 


91 


into continuum mechanics 

overcame the arbitrariness of the chosen rate in that it uniquely connects a particular (objective) strain rate 


37 Hooke’s law ; 1 lo] (cf. 141]) famously states that the strain in a deformation depends linearly on the occurring stress (“ut 
tensio, sic vis”). However, for finite deformations, different constitutive laws of elasticity can be obtained from this assumption, 
depending on the choice of a stress/strain pair. An idealized version of such a linear relation is given by ( |4.1l| ), i.e. by choosing 
the spatial Hencky strain tensor log V and the Kirchhoff stress tensor r. Since, however, Hooke speaks of extension versus force, 
the correct interpretati on of Hooke’s law is T Blot = 2 p (U — 1) + A tr {JJ — t) • 1, i.e. the case r = ^ in (TV2\ . 

38 Truesdell and Noll [203| p. 404] declared that “various such stress rates have been used in the literature. Despite claims and 
whole papers to the contrary, any advantage claimed for one such rate over another is pure illusion ”, and that “the properties of 
a material are independent of the choice of flux [i.e. of the chosen rate], which, like the choice of a [strain tensor], is absolutely 
immaterial ” 12031 p. 97]. 

39 For a shear test in Eulerian elasto-plasticity using the Zaremba-Jaumann rate ( |4.13| ), an unphysical artefact of oscillatory 
shear stress was observed, first in [122] . A similar oscillatory behavior was observed for hypoelasticity in [52] . 

40 Hill 108 used the terms conjugate and dual as synonyms. 
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to a stress tensor and, analogously, a stress rate to a strain tensor. The rational rule is that, when stress 
and strain tensors operate on configurations other than the reference configurations, the physically significant 
scalar products (S 2 , E\), {S 2 , Ef), (S 2 ,E\) and (S 2 , E\) (with the second Piola-Kirchhoff stress tensor S 2 and 
its work-conjugate Green strain tensor E\) must remain invariant, see 91 90 . 


4.2.2 Advantageous properties of the quadratic Hencky energy 

For modelling elastic material behavior there is no theoretical reason to prefer one strain tensor over another 
one, and the same is true for stress tensors. As discussed in Section |l.l[ stress and strain are immaterial^] 
Primary experimental data (forces, displacements) in material testing are sufficient to calculate any strain 
tensor and any stress tensor and to display any combination thereof in stress-strain curves, while only work- 
conjugate pairs are physically meaningful. 

However, for modelling finite-strain elasticity, the quadratic Hencky model 


Wh = H|dev„ log V\\ 2 
t = 2 fi dev„ log V + 


b 2 [ tr ( lo g v )f = t l ||dev„ log U|| 5 
k tr(log V) t , 


[tr(log U)] 2 


(4.15) 


exhibits a number of unique, favorable properties, including its functional simplicity and its dependency on 
only two material parameters pi and n that are determined in the infinitesimal strain regime and remain 
constant over the entire strain range. In view of the linear dependency of stress from logarithmic strain 
in (4.151, it is obvious that any nonlinearity in the stress-strain curves can only be captured in Hencky’s 


model by virtue of the nonlinearity in the strain tensor itself. There is a surprisingly large number of different 
materials, where Hencky’s elasticity relation provides a very good fit to experimental stress-strain data, which 
is true for different length scales and strain regimes. In the following we substantiate this claim with some 
examples. 

Nonlinear elasticity on macroscopic scales for a variety of materials. Anand j3,|4j has shown that the 
Hencky model is in good agreement with experiments on a wide class of materials, e.g. vulcanized natural 
rubber, for principal stretches between 0.7 and 1.3. More precisely, this refers to the characteristic that in 
tensile deformation the stiffness becomes increasingly smaller compared with the stiffness at zero strain, while 
for compressive deformation the stiffness becomes increasingly larger. 

Nonlinear elasticity in the very small strain regime. We mention in passing that a qualitatively similar 
dependency of material stiffness on the sign of the strain has been made much earlier in the regime of extremely 
small strains (10 -6 -10 -3 ). In Hartig’s law 89 from 1893 this dependency was expressed as = E° +ba, 


where E° is the elasticity modulus at zero stress and b < 0 is a dimensionless constant ]fjcf. the book of Bell 
19 and 126! in the context of linear elasticity with initial stress. Hartig also observed that the stress-stretch 
relation should have negative curvaturcf* 3 ] in the vicinity of the identity, as shown in Figure 14 

Crystalline elasticity on the nanoscale. Quite in contrast to the strictly stress-based continuum consti¬ 
tutive modelling, atomistic theories are based on a concept of interatomic forces. These forces are derived 
from potential^] V according to the potential relation f a = —d Xa V, which endows the model with a varia¬ 
tional structure. A further discussion of hybrid, atomistic-continuum coupling can be found in 60 . Thereby 


41 Cf. Truesdell 199, p. 145]: “It is important to realize that since each of the several material tensors [■ ■.] is an isotropic 
function of any one of the others, an exact description of strain in terms of any one is equivalent to a description in terms 
of any other” or Antman (o] p. 423]: “In place of C, any invertible tensor-valued function of C can be used as a measure of 
strain.” Rivlin jlTO] states that strain need never be defined at all, cf. |203[ p. 122]. 

42 The negative curvature (b < 0) was already suggested by Jacob Bernoulli in 1705 [2l] (cf. [ 20 ] p. 276]): “Homogeneous fibers 
of the same length and thickness, but loaded with different weights, neither lengthen nor shorten proportional to these weights; 
but the lengthening or the shortening caused by the small weight is less than the ratio that the first weight has to the second .” 

43 As Bell insists |l9| p. 155], a purely linear elastic response to finite strain, corresponding to zero curvature of the stress-strain 
curve at the identity 1, is never exhibited by any physical material: “The experiments of 280 years have demonstrated amply 
for every solid substance examined with sufficient care, that the [finite engineering] strain [U — t] resulting from small applied 
stress is not a linear function thereof.” 

44 For molecular dynamics (MD) simulations, a well-established level of sophistication is the modelling by potentials with 
environmental dependence (pair functionals like in the Embedded Atom Method (EAM) account for the energy cost to embed 
atomic nuclei into the electron gas of variable density) and angular dependence (like for Stillinger-Weber or Tersoff functionals). 
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the discreteness of matter at the nanoscale and the nonlocality of atomic interactions are inherently cap¬ 
tured. Here, atomistic stress is neither a constitutive agency nor does it enter a balance equation. Instead, it 
optionally can be calculated following the virial stress theorem 195 Chapter 8] to illustrate the state of the 
system. 


With their analyses in 53 and 54 , Dluzewski 


and coworkers aim to link the atomistic world to the 
macroscopic world of continuum mechanics. They 
search for the “best” strain measure with a view to¬ 
wards crystalline elasticity on the nanoscale. The 
authors consider the deformation of a crystal struc¬ 
ture and compare the atomistic and continuum ap¬ 
proaches. Atomistic calculations are made using 
the Stillinger-Weber potential. The stress-strain 
behaviour of the best-known anisotropic hyperelas¬ 
tic models are compared with the behaviour of the 
atomistic one in the uniaxial deformation test. The 
result is that the anisotropic energy based on the 
Hencky strain energy | (C. log U, log U), where C is 
the anisotropic elasticity tensor from linear elastic¬ 
ity, gives the best fit to atomistic simulations. More 
in detail, this best fit manifests itself in the observa¬ 
tion that for considerable compression (up to ~ 20%) 
the material stiffness is larger than the reference stiff¬ 
ness at zero strain, and for considerable tension (up 
to ss 20%) it is smaller than the zero-strain stiff¬ 
ness, again in good agreement with the atomistic result. This is also corroborated by comparing tabulated 
experimentally determined third order elastic constant^ 53 



Figure 14: The Biot stress T Blot corresponding to uniaxial 
stretches by factor A of incompressible materials fitted to 
experimental measurements by Jones and Treloar [115] . 
The curvature in A = 1 suggests negative third order con¬ 
stants (b < 0), which has also been postulated by Grioli 
[84| eq. (32)]. 


Elastic energy potentials based on logarithmic strain have also recently been motivated via molecular 
dynamics simulations 193 by Henann and Anand 194]. 


5 Applications and ongoing research 


5.1 The exponentiated Hencky energy 


As indicated in Section |TTT| and shown in Sections [2d] and [3] strain measures are closely connected to isotropic 
energy functions in nonlinear hyperelasticity: similarly to how the linear elastic energy may be obtained as 
the square of the Euclidean distance of Vw to so(n), the nonlinear quadratic Hencky strain energy is the 
squared Riemannian distance of Vy to SO(n). Fo r th e partial strain measures LOi so (F) = ||dev„ log VF T F\ 
and u m i(F) = |tr(log VF T F)\ defined in Theorem 


3.7 


the Hencky strain energy Wn can be expressed as 
Wu(F) = ^UF) + ^u J UF). (5.1) 


However, it is not at all obvious why this weighted squared sum should be viewed as the “canonical” energy 
associated with the geodesic strain measures: while it is reasonable to view the elastic energy as a quantity de¬ 
pending on some strain measure alone, the specific form of this dependence must not be determined by purely 
geometric deductions, but must take into account physical constraints as well as empirical observations]^] 

45 Third order elastic constants are corrections to the elasticity tensor in order to improve the response curves beyond the 
infinitesimal neighbourhood of the identity. They exist as tabulated values for many materials. Their numerical values depend 
on the choice of strain measure used which needs to be corrected. Dluzewski [53] shows that again the Hencky-strain energy 
4 (C. log U. log U) provides the best overall approx imat ion. 

46 G.W. Leibniz, in a letter to Jacob Bernoulli [l23| p. 572], stated as early as 1690 that “the [constitutive] relation between 
extension and stretching force should be determined by experiment ”, cf. 19, p. 10]. 
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Figure 16: The one-dimensional Hencky energy Wn compared to the exponentiated Hencky energy W e h and the corresponding 
Cauchy stresses <th, cr e H f° r very large uniaxial stretches A. Observe the non-convexity of Wh and the non-invertibility 
of <J H . 


For a large number of materials, the Hencky energy does indeed provide a very accurate model up to 
moderately large elastic deformations |3[|4|, i.e. up to stretches of about 40%, with only two constant material 
parameters which can be easily determined in the small strain range. For very large strain^] however, the 
subquadratic growth of the Hencky energy in tension is no longer in agreement with empirical measurements F] 
In a series of articles [154| 155, 153 80 , Neff et al. have therefore introduced the exponentiated Hencky energy 


W en {F) = + 

k 2k 


E e k II dev « lo S U \\ 2 _|_ — e k [tr(log U )] 2 

k 2k 


(5.2) 


with additional dimensionless material parameters 
f which for all values of k, k approx- 


k> \ and k > 


imates Hh for deformation gradients F sufficiently 
close to the identity 1, but shows a vastly different 
behaviour for ||F|| —» oo, cf. Figure 16 


The exponentiated Hencky energy has many ad¬ 
vantageous properties over the classical quadratic 
Hencky energy; for example, W e h is coercive on all 
Sobolev spaces W 1,p for 1 < p < oo, thus cavita¬ 
tion is excluded 12 143 . In the planar case n = 2, 
W e h is also polyconvex 155] 80 and thus Legendre- 
Hadamard-elliptic 10 , whereas the classical Hencky 


energy is not even LH-elliptic (rank-one convex) out¬ 
side a moderately large neighbourhood of 1 [36. [144 


(see also [113], where the loss of ellipticity for ener¬ 
gies of the form ||dev 3 log U\\P with hardening index 
0 < /3 < 1 are investigated). Therefore, many results 
guaranteeing the existence of energy-minimizing de¬ 
formations for a variety of boundary value problems 
can be applied directly to W e h for n = 2. 



Figure 15: The equation of state (EOS), i.e. the trace of the 
Cauchy stress corresponding to a purely volumetric de¬ 
formation (cf. 168] ), for the quadratic and the exponen¬ 
tiated Hencky model (with parameter k = 4). 


47 The elastic range of numerous materials, including vulcanized rubber or skin and other soft tissues, lies well above stretches 
of 40%. 

48 While the behaviour of elasticity models for extremely large strains might not seem important due to physical restraints and 
intermingling plasticity effects outside a narrow range of perfect elasticity, it is nevertheless important to formulate an idealized 
law of elasticity over the whole range of deformations; cf. Hencky [99] p. 215] (as translated in [146] p.2]): “It is not important 
that such an idealized elastic [behaviour] does not actually exist and our ideally elastic material must therefore remain an ideal. 
Like so many mathematical and geometric concepts, it is a useful ideal, because once its deducible properties are known it can 
be used as a comparative rule for assessing the actual elastic behaviour of physical bodies 
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127 


107 


Furthermore, W e h satisfies a number of constitutive inequalities 154 such as the Baker-Ericksen inequality 
, the pressure-compression inequality and the tension-extension inequality as well as Hill’s inequalitjj^] 


160 161 , which is equivalent to the convexity of the elastic energy with respect to the logarithmic strain 


tensor 192 


Moreover, for W e h, the Cauchy-stress-stretch relation V <r e H(V) is invertible (a property hitherto 

unknown for other hyperelastic f ormu lations) and pure Cauchy shear stress corr esponds t o pure shear strain, 


as is the case in linear elasticity 


154 


The physical meaning of Poisson’s ratio 


169 


79 


V 2(3K+£t) a ^ S0 

similar to the linear case; for example, v = 1 directly corresponds to incompressibility of the material and 
v = 0 implies that no lateral extension or contraction occurs in uniaxial tensions tests. 


5.2 Related geodesic distances 

The logarithmic distance measures obtained in Theorems |3.3| ancl |3.7| show a strong similarity to other geodesic 
distance measures on Lie groups. For example, consider the special orthogonal group SO(n) endowed with 
the canonical bi-invariant Riemannian metrical 

9q(X,Y) = (Q t X,Q t Y) = (X,Y) 

for Q G SO(n) and A', Y G Tq SO(n) = Q ■ so(n). Then the geodesic exponential at t G SO(n) is given by the 
matrix exponential on the Lie algebra so(n), i.e. all geodesic curves are one-parameter groups of the form 

7 (t) = Q ■ exp [tA) 

with Q G SO(ra) and A G so(n) (cf. 136|). It is easy to show that the geodesic distance between Q,R G SO (n) 
with respect to this metric is given by 

distg eoc i ; SO(n) (Q;-^) ||log((5 ^)ll ^ 

where || . || is the Frobenius matrix norm and log: SO(n) —> so (n) denotes the principal matrix logarithm on 
SO(n), which is uniquely defined by the equality exp(logQ) = Q and the requirement Ai(logQ) G (—7r, 7r] for 
all Q G SO(n) and all eigenvalues Ai(logQ). 

This result can be extended to the geodesic distance on the conformal special orthogonal group CSO(ro) 
consisting of all angle-preserving linear mappings: 

CSO(n) := {c • Q \ c > 0, Q G SO(n)} , 


where the bi-invariant metric 5csO(n) is given by the canonical inner product: 

g° so{n) (X,Y) = (A~ 1 X,A~ 1 Y ). 

Then 

distgeod.csowtc- Q,d-R) = ||log(<2 i ‘i?)f + 


(5.3) 


lr -(r 


n l 


where log again denotes the principal matrix logarithm on SO(n). Note that the punctured complex plane 
C \ {0} can be identified with CSO(2) via the mapping 


z = a + ib i ^ Z G CSO(2) = 


a 

-6 


a 2 + b 2 ± 0 


161 


can be stated more generally as ( “.l [r] — m [r D — D r], D) > 0 in the hypoelastic formulation, where 


49 Hill’s inequality 

d ° is the Zaremba-Jaumann stress rate and r is the Kirchhoff stress tensor. For m = 0, as Silhavy explains, “Hill’s 


d t 


inequalities [. ..] require the convexity of [the strain ener gy W] in [terms of the strain tensor log V] . . . This does not seem to 
contradict any theoretical or experimental evidence ” 1931 p. 309]. 

50 Note that p c • is the restriction of our left-GL(n)- invariant, right-0(n)-invariant metric g (as defined in Section |3.l| ) to 
SO(n). 
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5.3 Outlook 


While first applications of the exponentiated Hencky energy, which is based on the partial strain measures 
u! iso, w vo i introduced here, show promising results, including an accurate modelling of so-called tire-derived 
material 140 , a more thorough fitting of the new parameter set to experimental data is necessary in order 
to assess the range of applicability of W e yi towards elastic materials like vulcanized rubber. A different 
formulation in terms of the partial strain measures Wj so and w vo i, i.e. an energy function of the form 


W(F) = ^(uj iso (F),u} vo i(F)) = 'F(||dev 3 log£/||, |tr(log[/)|) 


(5.4) 


with d': [0, oo ) 2 —► [0, oo), might even prove to be polyconvex in the three-dimensional case. The main open 
problem of finding a polyconvex (or rank-one convex) isochoric energy function 51 ] A i —> \k(||dev 3 log t/||) has 


also been considered by Sendova and Walton 189 . Note that while every isotropic elastic energy W can be 


expressed as W(F) = h(Ki, K 2 , A 3 ) with Criscione’s invariant 


I< 1 = tr(log U ), K 2 = ||dev 3 log U\\ and 


52 

46 

45 

51 

208 


K, 


/ dev 3 log U \ 
6 \ ||dev 3 log U\\) 


(5.5) 


not every elastic energy has a representation of the form (5.4); for example, (5.4) implies the tension- 
compression symmetr^\W (F) = W(F ~ 1 ), which is not necessarily satisfied by energy functions in general p^] 
In terms of the Shield transformation 55 191 39 


W*{F) := detF • W{F ~ l ), 

the tension-compression symmetry amounts to the requirement dc 5 F W* (F) = W ( F ) or, for incompressible 
materials, W*(F) = W(F). Moreover, under the assumption of incompressibility, the symmetry can be 
immediately extended to arbitrary deformations <p: fl — > <p(0) and : V 5 )^) ~^ fJ: if det V<p = 1, we can 
apply the substitution rule to find 



W(V(y> 1 )(x))da; 


l )(^{x))) • |det Vy>(a;)| da; 


W(V<p( x) 1 )dx 


W{\/ip{x)) dec 


if W(F 1 ) = W(F) for all F £ SL(n), thus the total energies of the deformations ip,<p 1 are equal, cf. Figure 

m 

Since the function 

P |_). e l|de5V 2 log c/|| 2 _ g dist geod,SL(2) ( det ^V2 ’ SO ( 2 )) 

51 Ideally, the function ^ should also satisfy additional requirements, such as monotonicity, convexity and exponential growth. 

52 The invariants K\ and K% = tr ((dev 3 log U) 2 ) as well as K% = tr ((dev 3 log(7) 3 ) had already been discussed exhaustively 
by H. Richter in a 1949 ZAMM article [176| §4], while K\ and Ki have also been considered by A.I. Lurie |125| p. 189]. Criscione 
has shown that the invariants given in ( |5.5| ) enjoy a favourable orthogonality condition which is useful when determining material 
parameters. 

53 The tension-compression symmetry is often expressed as r(V~ 1 ) = — t(V), where t(V) is the Kirchhoff stress tensor 
corresponding to the left Biot stretch V . This condition, which is the natural nonlinear counterpart of the equality cr(— e) = — cr(e) 
in linear elasticity, is eq uival ent to the condition W(F~ 1 ) = W(F) for hyperelastic constitutive models. 

54 Truesdell and Noll [203[ p. 174] argue that “.. . there is no foundation for the widespread belief that according to the theory 
of elasticity, pressure and tension have equal but opposite effects ”. Examples for isotropic energy functions which do not satisfy 
this symmetry condition in general but only in the incompressible case can be found in [92] . For an idealized isotropic elastic 
material, however, the tension-compression symmetry is a plausible requirement (with an obvious additive counterpart in linear 
elasticity), especially for incompressible bodies. 

55 Further properties of the Shield transformation can be found in 193| p.288]; for example, it preserves the poly convexity, 
quasiconvexity and rank-one convexity of the original energy. 
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Figure 17: The tension-compression symmetry for incompressible materials: if detV</? = 1 and W(F 1 ) = W(F) for all 
F £ SL(n), then f n W(Vip(x)) dx = W(V(vJ- 1 )(a:)) dx. 


in planar elasticity is polyconvex 155 80 , it stands to reason that a similar formulation in the three- 


dimensional case might prove to be polyconvex as well. A first step towards finding such an energy is to 
identify where the function W with 


W(F) = e ll dev 3 l° g £d| 2 = g dlSt geod, SL(3) ( det pl/3 ’ SO ( 3 )) 


(5.6) 

which is not rank-one convex [154] , loses its ellipticity properties. For that purpose, it may be useful to 
consider the quasiconvex hull of W. There already are a number of promising results for similar energy 
functions; for example, the quasiconvex hull of the mapping 

F | -t distE U ciidT) SO(2)) = \\U — 1 || 2 


can be explicitly computed |194j 56 IHZl , and the quasiconvex hull of the similar Saint-Venant-Kirchhoff energy 

T/T7). -r-, r ( 7?'l LL. 11 _ 'll I I 2 I ^ fTr- ^ 'll M 2 V» o o rmrnri t T TAy - F oni) T? onnlf 101 T-TVvv +- V>/-» m oVI rro 


Wsvk(F) = 4 \\C — H || 2 + | [tr(C — I )] 2 has been given by Le Dret and Raoult 121 . For the mappings 


F i-a dist| uclid (F, SO(3)) or F>-> dist geod (F, SO(ra)) 
with n > 2 , however, no explicit representation of the quasiconvex hull is yet known, although it has been 


shown that both expressions are not rank-one convex 24 


It might also be of interest to calcul ate the geodesic distance dist geo d(A, B ) for a larger class of matrices 


A,B£ GL + (n)j^| although Theorem 3.3 allows us to explicitly compute the distance dist geo d(l, P) for 
P £ Sym + (n) and local results are available for certain special cases m. it is an open question whether 
there is a general formula for the distance dist geodi G L+(n) (Qi R) between arbitrary rotations R, Q £ SO(n) for 
all parameters fi, k > 0. Since restricting our left-GL(n)-invariant, right-0(n)-invariant metric on GL(ro) 
to SO(n) yields a multiple of the canonical bi-SO(n)-invariant metric on SO(n), we can compute 

dist leod,Gh+(n)(Qi R ) = Me ' dist g eod , SO(n) (Q, i?) = Me ||fog(Q T 1 ?) || 2 

if for all Q,R £ SO(n) a shortest geodesic in GL + (n) connecting Q and R is already contained within SO(n), 
cf. Figure [Ts} However, whether this is the case depends on the chosen parameters fj,, q c ; a general closed-form 
solution for dist geodj gl+O) on SO(n) is therefore not yet known 


128 


Moreover, it is not known whether our result can be generalized to anisotropic Riemannian metrics, 
i.e. if the geodesic distance to SO(n) can be explicitly computed for a larger class of left-GL(n)-invariant 


56 An improved understanding of the geometric structure of mechanical problems could, for example, help to develop new 
discretization methods 184 85], 
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Figure 18: If SO(n) contains a length minimizing geodesic connecting Q,R E SO(n) with respect to our left-GL(n)-invariant, 
right-0(n)-invariant metric g on GL(n), then the GL - *" (n)-geodesic distance between Q and R is equal to the well-known 
SO(n)-geodesic distance /i c ||log(Q T i?)|| 2 . 


Riemannian metrics which are not necessarily right-O(n)-invariant. A result in this direction would have 
immediate impact on the modelling of finite strain anisotropic elasticity 14 188 187 . The difficulties with 


such an extension are twofold: one needs a representation formula for Riemannian metrics which are right- 
invariant under a given symmetry subgroup of O(n), as well as an understanding of the corresponding geodesic 


curves. 


6 Conclusion 

We have shown that the squared geodesic distance of the (finite) deformation gradient F £ GL + (n) to the 
special orthogonal group SO(n) is the quadratic isotropic Hencky strain energy: 

distg eod (F, SO(n)) = /z ||dev n log U\\ 2 + ^ [tr(log U )] 2 , 

if the general linear group is endowed with the left-GL(ro)-invariant, right-0(?z)-invariant Riemannian metric 
g A (X,Y) = {A- 1 X,A~ 1 Y) li>lictK , where 

(X, Y)= n (dev„ sym X, dev n sym Y) + /z c (skew X, skew Y) + § tr(X) tr(Y) 

with (X, Y) = tr (X T Y). Furthermore, the (partial) logarithmic strain measures 

w iso = ||dev n log U\\ = ||dev ra log VF T F\\ and w vo i = |tr(log C/)| = |tr(log VF T F)\ 

have been characterized as the geodesic distance of F to the special orthogonal group SO(n) and the identity 
tensor 1, respectively: 

w iso = |dev n log U\\ = distgeod, SL(n) ^ pVn ’ > 

Wvoi = |tr(log U)\ = yjn- dist geodi K+-! ((detF) 1/n • 1, l) , 

where the geodesic distances on SL(n) and R + -l are induced by the canonical left invariant metric g A (X, Y) = 
(A~ l X,A~ 1 Y). 

We thereby show that the two quantities Wi so = ||dev ra log U\\ and w vo i = |tr(log?7)| are purely geometric 
properties of the deformation gradient F, similar to the invariants ||dev n e|| and |tr(e)| of the infinitesimal 
strain tensor e in the linearized setting. 

While there have been prior attempts to deductively motivate the use of logarithmic strain in nonlinear 
elasticity theory, these attempts have usually focussed on the logarithmic Hencky strain tensor Eq = log U 
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(or Eg = log V") and its status as the “natural” material (or spatial) strain tensor in isotropic elasticity. We 
discussed, for example, a well-known characterization of logH in the hypoelastic context: if the strain rate 
^ is objective as well as corotational , and if 

d D 

^ [E] = D := sym(FF" 1 ) 

for some strain tensor E, then ^ = 4*° s must be the logarithmic rate and E = Eq = \ogV must be the 
spatial Hencky strain tensor. 

However, as discussed in Section |1.1[ all strain tensors are interchangeable: the choice of a specific strain 
tensor in which a constitutive law is to be expressed is not a restriction on the available constitutive relations. 
Such an approach can therefore not be applied to deduce necessary conditions or a priori properties of 
constitutive laws. 

Our deductive approach, on the other hand, directly motivates the use of the strain measures Wi so and w vo i 
from purely differential geometric observations. As we have indicated, the requirement that a constitutive 
law depends only on Wj so and w V oi has direct implications; for example, the tension-compression symmetry 
W(F) = W{F~ 1 ) is satisfied by every hyperelastic potential W which can be expressed in terms of Wi so and 
Wvoi alone. 

Moreover, as demonstrated in Section [4} similar approaches oftentimes presuppose the role of the positive 
definite factor U = V F T F as the sole measure of the deformation, whereas this independence from the 
orthogonal polar factor is obtained deductively in our approach (cf. Table [l]). 



Measure of deformation deduced 

Measure of deformation postulated 

linear 

dist| uc i id ^^ C!K (VM,so(n)) 

= g ||dev n symVw || 2 + ^ [tr(sym Vtt )] 2 

dist Euclid, 0) 

= g ||dev n e 2 + [tr(e )] 2 

geometrically 

nonlinear 

dist| uclid (T, SO(n)) = g \\VWf - 1|| 2 

disthuclid,Sym(n)(G, t) = g \\U — t\\~ 

geometrically 

nonlinear 

(weightedjp 7 ] 

not well defined 

dist Euc ii d , Sym(n),/i,(t(^i ^-) 

= \x ||dev n (£/ — 1 )|| 2 + ^ [tr(C7 — l )] 2 

geodesic 

distgeod,, SO(?l)) 

= g ||dev n log(VF T F)\\ 2 + ^ [tr(log VF T F)] 2 

dist geodi 

= g ||dev n log G || 2 + | [tr(log U)] 2 

log-Euclidean 

not well defined 

dist Log-Euclid,/^,/c (^j l) 

distEuciid, Sym(n),/x,«(^S 0) 

= g ||dev n log U\\ 2 + | [tr(log U)} 2 


Table 1: Different approaches towards the motivation of different strain tensors and strain measures. 


Note also that the specific distance measure dist geo d on GL + (n) used here is not chosen arbitrarily: the 
requirements of left-GL(n)-invariance and right-0(n)-invariance, which have been motivated by mechanical 
considerations, uniquely determine g up to the three parameters g, g Cl k > 0. This uniqueness property further 
emphasizes the generality of our results, which yet again strongly suggest that Hencky’s constitutive law 


57 Observe that |r]ev,, (f/ — 1)|| 2 does not measure the isochoric (distortional) part 


F 

(detF) 1 /’ 


of F. 
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should be considered the idealized nonlinear model of elasticity for very small strains outside the infinitesimal 
range. 


Acknowledgements 


The second author acknowledges support by the Deutsche Forschungsgemeinschaft (DFG) through a Heisen¬ 
berg fellowship under grant El 453/2-1. 

We are grateful to Prof. Alexander Mielke (Weierstrafi-Institut, Berlin) for pertinent discussions on 
geodesics in GL(ra); the first parametrization of geodesic curves on SL(n) known to us is due to him 134 . We 


also thank Prof. Robert Bryant (Duke University) for his helpful remarks regarding geodesics on Lie groups 
and invariances of inner products on 0 l(n), as well as a number of friends who helped us with the draft. 

We also thank Dr. Andreas Fischle (Technische Universitat Dresden) who, during long discussions on 
continuum mechanics and differential geometry, inspired many of the ideas laid out in this paper. 

The first author had the great honour of presenting the main ideas of this paper to Richard Toupin on the 
occasion of the Canadian Conference on Nonlinear Solid Mechanics 2013 in the mini-symposium organized 
by Francesco dcll’Isola and David J. Steigmann, which was dedicated to Toupin. 


Conflict of Interest 

The authors declare that they have no conflict of interest. 

References 

[1] A. H. Al-Mohy, N. J. Higham, and S. D. Relton. “Computing the Frechet derivative of the matrix logarithm and 
estimating the condition number”. SIAM Journal on Scientific Computing 35.4 (2013), pp. C394-C410. 

[2] E. Almansi. “Sulle deformazioni finite dei solidi elastici isotropi”. Rendiconti della Reale Accademia dei Lincei, Classe 
di scienze fisiche, matematiche e naturali 20 (1911). 

[3] L. Anand. “On H. Hencky’s approximate strain energy function for moderate deformations”. Journal of Applied Me¬ 
chanics 46 (1979), pp. 78-82. 

[4] L. Anand. “Moderate deformations in extension-torsion of incompressible isotropic elastic materials”. Journal of the 
Mechanics and Physics of Solids 34 (1986), pp. 293-304. 

[5] E. Andruchow, G. Larotonda, L. Recht, and A. Varela. “The left invariant metric in the general linear group”. Journal 
of Geometry and Physics 86.0 (2014), pp. 241 -257. 

[6] S. S. Antman. Nonlinear problems of elasticity. Vol. 107. Applied Mathematical Sciences. New York: Springer, 2005. 

[7] V. Arsigny, O. Commowick, N. Ayache, and X. Pennec. “A fast and log-Euclidean polyaffine framework for locally linear 
registration”. Journal of Mathematical Imaging and Vision 33.2 (2009), pp. 222-238. 

[8] V. Arsigny, P. Fillard, X. Pennec, and N. Ayache. “Fast and simple calculus on tensors in the log-Euclidean framework”. 
Medical Image Computing and Computer-Assisted Intervention-MICCAI 2005. Springer, 2005, pp. 115-122. 

[9] V. Arsigny, P. Fillard, X. Pennec, and N. Ayache. “Geometric means in a novel vector space structure on symmetric 
positive-definite matrices”. SIAM Journal on Matrix Analysis and Applications 29.1 (2007), pp. 328-347. 

[10] J. M. Ball. “Convexity conditions and existence theorems in nonlinear elasticity”. Archive for Rational Mechanics and 
Analysis 63.4 (1976), pp. 337-403. 

[11] J. M. Ball. “Constitutive inequalities and existence theorems in nonlinear elastostatics”. Nonlinear Analysis and Me¬ 
chanics: Heriot-Watt Symposium. Ed. by R. Knops. Vol. 1. Pitman Publishing Ltd. Boston. 1977, pp. 187-241. 

[12] J. M. Ball. “Discontinuous equilibrium solutions and cavitation in nonlinear elasticity”. Philosophical Transactions of 
the Royal Society of London. Series A, Mathematical and Physical Sciences 306.1496 (1982), pp. 557-611. 

[13] J. M. Ball. “Some open problems in elasticity”. Geometry, Mechanics, and Dynamics. Ed. by P. Newton, P. Holmes, 
and A. Weinstein. Springer, 2002, pp. 3-59. 

[14] D. Balzani, P. Neff, J. Schroder, and G. A. Holzapfel. “A polyconvex framework for soft biological tissues. Adjustment 
to experimental data”. International Journal of Solids and Structures 43.20 (2006), pp. 6052-6070. DOI: 10.1016/j . 
ij solstr.2005.07.048 


38 



[15] R. C. Batra. “Linear constitutive relations in isotropic finite elasticity”. Journal of Elasticity 51.3 (1998), pp. 243-245. 

[16] R. C. Batra. “Comparison of results from four linear constitutive relations in isotropic finite elasticity”. International 
Journal of Non-Linear Mechanics 36.3 (2001), pp. 421-432. 

[17] Z. P. Bazant. “Easy-to-compute tensors with symmetric inverse approximating Hencky finite strain and its rate”. Journal 
of Engineering Materials and Technology 120.2 (1998), pp. 131-136. 

[18] G. F. Becker. “The Finite Elastic Stress-Strain Function”. American Journal of Science 46 (1893). newly typeset version 
available at https://www.uni-due.de/imperia/md/content/matheinatik/ag_neff/becker_latex_newl893.pdf, pp. 337- 
356. 

[19] J. Bell and C. Truesdell. Mechanics of Solids: Volume 1: The Experimental Foundations of Solid Mechanics. Handbuch 
der Physik. Springer, 1973. 

[20] E. Benvenuto. An Introduction to the History of Structural Mechanics. Part I: Statics and Resistance of Solids. Springer, 
1991. 

[21] J. Bernoulli. “Veritable hypothese de la resistance des solides, avec la demonstration de la courbure des corps qui font 
ressort”. Memoires de I’Academie des Sciences (1705). 

[22] D. S. Bernstein. Matrix Mathematics: Theory, Facts, and Formulas (Second Edition). Princeton reference. Princeton 
University Press, 2009. 

[23] A. Bertram. Elasticity and Plasticity of Large Deformations. Springer, 2008. 

[24] A. Bertram, T. Bohlke, and M. Silhavy. “On the rank 1 convexity of stored energy functions of physically linear stress- 
strain relations”. Journal of Elasticity 86.3 (2007), pp. 235-243. 

[25] R. Bhatia. Positive definite matrices. Princeton University Press, 2009. 

[261 R. Bhatia and J. Holbrook. “Riemannian geometry and matrix geometric means”. Linear Alqebra and its Applications 
413.2 (2006), pp. 594-618. 

[27] D. Bigoni. Nonlinear Solid Mechanics: Bifurcation Theory and Material Instability. Cambridge University Press, 2012. 

[28] M. A. Biot. “Non-linear theory of elasticity and the linearized case for a body under initial stress”. The London, 
Edinburgh, and Dublin Philosophical Magazine and Journal of Science 27.183 (1939), pp. 468-489. 

[29] M. Birsan, P. Neff, and J. Lankeit. “Sum of squared logarithms - an inequality relating positive definite matrices and their 
matrix logarithm”. Journal of Inequalities and Applications 2013.1 (2013), pp. 1-16. DOI: 10.1186/1029-242X-2013-168 

[30] L. Borisov, P. Neff, S. Sra, and C. Thiel. “The sum of squared logarithms inequality in arbitrary dimensions”, to appear 
in Linear Algebra and its Applications (2015). available at arXiv: 1508.04039. 

[31] C. Bouby, D. Fortune, W. Pietraszkiewicz, and C. Vallee. “Direct determination of the rotation in the polar decom¬ 
position of the deformation gradient by maximizing a Rayleigh quotient”. Zeitschrift fur Angewandte Mathematik und 
Mechanik 85.3 (2005), pp. 155-162. 

[32] R. Brannon. Define your strain! available at www.mech.utah.edu/~brannon/public/strain.pdf 

[33] O. T. Bruhns. “Some remarks on the history of plasticity - Heinrich Hencky, a pioneer of the early years”. The History 
of Theoretical, Material and Computational Mechanics - Mathematics Meets Mechanics and Engineering. Ed. by E. 
Stein. Springer, 2014, pp. 133-152. 

[341 O. T. Bruhns. “The Prandtl-Reuss equations revisited”. Zeitschrift fur Anqewandte Mathematik und Mechanik 94.3 
(2014), pp. 187-202. 

[35] O. T. Bruhns, A. Meyers, and H. Xiao. “On non-corotational rates of Oldroyd’s type and relevant issues in rate con¬ 
stitutive formulations”. Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences 460.2043 
(2004), pp. 909-928. DOI: 10.1098/rspa.2003.1184 

[36] O. T. Bruhns, H. Xiao, and A. Mayers. “Constitutive inequalities for an isotropic elastic strain energy function based on 
Hencky’s logarithmic strain tensor”. Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences 
457 (2001), pp. 2207-2226. 

[37] O. T. Bruhns. “The multiplicative decomposition of the deformation gradient in plasticity—origin and limitations”. 
From Creep Damage Mechanics to Homogenization Methods - A Liber Amicorum to celebrate the birthday of Nobutada 
Ohno. Ed. by H. Altenbach, T. Matsuda, and D. Okumura. Vol. 64. Springer International Publishing, 2015, pp. 37-66. 

[38] M. M. Carroll. “Must elastic materials be hyperelastic?” Mathematics and Mechanics of Solids 14.4 (2009), pp. 369-376. 

[39] M. M. Carroll and F. J. Rooney. “Implications of Shield’s inverse deformation theorem for compressible finite elasticity”. 
Zeitschrift fur angewandte Mathematik und Physik 56.6 (2005), pp. 1048-1060. 

[40] A. L. Cauchy. Cours d’analyse de VEcole royale poly technique: I. Analyse algebrique. available at https://archive. 
org/details/coursdanalysedeOOcaucgoog Paris: Impr. royale Debure freres, 1821. 

[41] A. L. Cauchy. “Sur la condensation et la dilatation des corps solides”. Exercices de Mathematiques. Vol. 2. Chez de 
Bure freres, 1827, pp. 60-69. 


39 



[42] A. L. Cauchy. “Memoire sur les dilatations, les condensations et les rotations produits par un changement de forme dans 
un systeme de points materiels”. CEuvres completes d’Augustin Cauchy. Vol. XII. available at http://gallica.bnf.fr/ 
ark:/12148/bpt6k90204r/f346 Gauthier-Villars, 1841. 

[43] I. Chao, U. Pinkall, P. Sanan, and P. Schroder. “A simple geometric model for elastic deformations”. ACM Transactions 
on Graphics. Vol. 29. 4. ACM. 2010, 38:1-38:6. 

[44] P. G. Ciarlet. Three-Dimensional Elasticity. Studies in Mathematics and its Applications 1. Elsevier Science, 1988. 

[45] J. C. Criscione. “Direct tensor expression for natural strain and fast, accurate approximation”. Computers and Structures 
80.25 (2002), pp. 1895-1905. 

[46] J. C. Criscione, J. D. Humphrey, A. S. Douglas, and W. C. Hunter. “An invariant basis for natural strain which yields 
orthogonal stress response terms in isotropic hyperelasticity”. Journal of the Mechanics and Physics of Solids 48.12 
(2000), pp. 2445-2465. 

[47] A. Curnier and L. Rakotomanana. “Generalized strain and stress measures: critical survey and new results”. Engineering 
Transactions, Polish Academy of Sciences 39.3-4 (1991), pp. 461-538. 

[48] F. M. Dannan, P. Neff, and C. Thiel. “On the sum of squared logarithms inequality and related inequalities”, to appear 
in JMI: Journal of Mathematical Inequalities (2014). availabelt at arXiv:1411.1290. 

[49] H. Darijani and R. Naghdabadi. “Constitutive modeling of solids at finite deformation using a second-order stress-strain 
relation”. International Journal of Engineering Science 48.2 (2010), pp. 223 -236. 

[50] C. De Boor. “A naive proof of the representation theorem for isotropic, linear asymmetric stress-strain relations”. Journal 
of Elasticity 15.2 (1985), pp. 225-227. ISSN: 0374-3535. DOI: 10.1007/BF00041995 

[51] J. Diani and P. Gilormini. “Combining the logarithmic strain and the full-network model for a better understanding 
of the hyperelastic behavior of rubber-like materials”. Journal of the Mechanics and Physics of Solids 53.11 (2005), 
pp. 2579-2596. 

[52] J. Dienes. “On the analysis of rotation and stress rate in deforming bodies”. Acta Mechanica. 32 (1979), pp. 217-232. 

[53] P. Dluzewski. “Anisotropic hyperelasticity based upon general strain measures”. Journal of Elasticity 60 (2000), pp. 119- 
129. 

[54] P. Dluzewski and P. Traczykowski. “Numerical simulation of atomic positions in quantum dot by means of molecular 
statics.” Archives of Mechanics 55 (2003), pp. 501—514. 

[55] A. Dolcetti and D. Pertici. “Some differential properties of GL n (R) with the trace metric”. Rivista di Matematica della 
Universita di Parma 6 (2015), pp. 267-286. 

[56] G. Dolzmann. “Regularity of minimizers in nonlinear elasticity - the case of a one-well problem in nonlinear elasticity”. 
Technische Mechanik 32 (2012), pp. 189-194. 

[57] G. Dolzmann, J. Kristensen, and K. Zhang. “BMO and uniform estimates for multi-well problems”. Manuscripta Math- 
ematica 140.1-2 (2013), pp. 83-114. 

[58] T. C. Doyle and J. L. Ericksen. “Nonlinear elasticity”. Advances in Applied Mechanics 4 (1956), pp. 53-115. 

[59] B. Eidel, P. Neff, and R. J. Martin. “Tractatus mathematicus-mechanicus modis deformationibus sub aspecto geometriae 
differentialis”. in preparation (2016). 

[60] B. Eidel and A. Stukowski. “A variational formulation of the quasicontinuum method based on energy sampling in 
clusters”. Journal of the Mechanics and Physics of Solids 57.1 (2009), pp. 87-108. 

[61] M. Epstein. The Geometrical Language of Continuum Mechanics. Cambridge University Press, 2010. 

[62] L. Euler. Methodus inveniendi lineas curvas maximi minimive proprietate gaudentes (appendix, de curvis elasticis). 
Lausannae & Genevae, 1774. 

[63] S. Federico. “Some Remarks on Metric and Deformation”. Mathematics and Mechanics of Solids 20 (2015), pp. 522-539. 

[64] Z. Fiala. “Time derivative obtained by applying the Riemannian manifold of Riemannian metrics to kinematics of 
continua”. Comptes Rendus Mecanique 332.2 (2004), pp. 97-102. 

[65] Z. Fiala. “Geometrical setting of solid mechanics”. Annals of Physics 326.8 (2011), pp. 1983-1997. 

[66] Z. Fiala. “Evolution equation of Lie-type for finite deformations, time-discrete integration, and incremental methods”. 
Acta Mechanica 226.1 (2014), pp. 17-35. 

[67] Z. Fiala. “Discussion of ”On the interpretation of the logarithmic strain tensor in an arbitrary system of representation“ 
by M. Latorre and F.J. Montans”. International Journal of Solids and Structures 56-57.0 (2015), pp. 290 -291. 

[68] Z. Fiala. “Geometry of finite deformations and time-incremental analysis”. International Journal of Non-Linear Me¬ 
chanics 81 (2016), pp. 230-244. 


40 


[69] J. Finger. “Das Potential der inneren Krafte und die Beziehungen zwischen den Deformationen und den Spannungen 
in elastisch isotropen Korpern bei Beriicksichtigung von Gliedern, die beziiglich der Deformationselemente von dritter, 
beziehungsweise zweiter Ordnung sind”. Sitzungsberichte der Akademie der Wissenschaften in Wien 44 (1894). ISSN: 
0378-8644. 

[70] A. Fischle and P. Neff. “The geometrically nonlinear Cosserat micropolar shear-stretch energy. Part I: A general parame¬ 
ter reduction formula and energy-minimizing microrotations in 2D”. to appear in Zeitschrift fur Angewandte Mathematik 
und Mechanik (2015). available at arXiv: 1507.05480. 

[71] A. Fischle and P. Neff. “The geometrically nonlinear Cosserat micropolar shear-stretch energy. Part II: Non-classical 
energy-minimizing microrotations in 3D and their computational validation”, submitted (2015). available at arXiv: 1509.06236. 

[72] J. E. Fitzgerald. “A tensorial Hencky measure of strain and strain rate for finite deformations”. Journal of Applied 
Physics 51.10 (1980), pp. 5111-5115. 

[73] P. J. Flory. “Thermodynamic relations for high elastic materials”. Transactions of the Faraday Society 57 (1961), 
pp. 829-838. 

[74] R. L. Fosdick and A. S. Wineman. “On general measures of deformation”. Acta Mechanica 6.4 (1968), pp. 275—295. 

[75] A. D. Freed. “Natural strain”. Journal of Engineering Materials and Technology 117.4 (1995), pp. 379-385. 

[76] A. D. Freed. “Hencky strain and logarithmic rates in Lagrangian analysis”. International Journal of Engineering Science 
81 (2014), pp. 135-145. 

[77] A. D. Freed. Soft Solids - A Primer to the Theoretical Mechanics of Materials. Birkhauser Basel, 2014. 

[78] G. Friesecke, R. D. James, and S. Muller. “A theorem on geometric rigidity and the derivation of nonlinear plate theory 
from three-dimensional elasticity”. Communications on Pure and Applied Mathematics 55.11 (2002), pp. 1461-1506. 

[79] H. Gercek. “Poisson’s ratio values for rocks”. International Journal of Rock Mechanics and Mining Sciences 44.1 (2007), 
pp. 1-13. 

[80] I.-D. Ghiba, P. Neff, and M. Silhavy. “The exponentiated Hencky-logarithmic strain energy. Improvement of planar 
polyconvexity”. International Journal of Non-Linear Mechanics 71 (2015), pp. 48-51. DOI: 10.1016/j . ijnonlinmec . 
2015.01.009 

[81] G. Green. “On the propagation of light in crystallized media”. Transactions of the Cambridge Philosophical Society 7 
(1841), p. 121. 

[82] G. Grioli. “Una proprieta di minimo nella cinematica delle deformazioni finite.” Bollettino deWUnione Matematica 
Italiana 2 (1940), pp. 252-255. 

[83] G. Grioli. Mathematical Theory of Elastic Equilibrium (recent results). Vol. 7. Ergebnisse der angewandten Mathematik. 
Springer, 1962. 

[84] G. Grioli. “On the thermodynamic potential for continuums with reversible transformations - some possible types”. 
Meccanica 1.1-2 (1966), pp. 15-20. 

[85] P. Grohs, H. Hardering, and O. Sander. “Optimal a priori discretization error bounds for geodesic finite elements”. 
Foundations of Computational Mathematics (2013), pp. 1—55. 

[86] M. E. Gurtin and K. Spear. “On the relationship between the logarithmic strain rate and the stretching tensor”. 
International Journal of Solids and Structures 19.5 (1983), pp. 437-444. 

[87] K. Hackl, A. Mielke, and D. Mittenhuber. “Dissipation distances in multiplicative elastoplasticity”. Analysis and Sim¬ 
ulation of Multifield Problems. Ed. by W. L. Wendland and M. Efendiev. Springer, 2003, pp. 87-100. 

[88] M. Hanin and M. Reiner. “On isotropic tensor-functions and the measure of deformation”. Zeitschrift fiir angewandte 
Mathematik und Physik 7.5 (1956), pp. 377-393. 

[89] E. Hartig. “Der Elastizitatsmodul des geraden Stabes als Funktion der spezifischen Beanspruchung”. Der Civilingenieur 
39 (1893). available at www.uni-due.de/imperia/md/content/mathernatik/ag_neff/hartig_elastizitaetsmodul.pdf, 
pp. 113-138. 

[90] P. Haupt. Continuum Mechanics and Theory of Materials. Heidelberg: Springer, 2000. 

[91] P. Haupt and C. Tsakmakis. “On the application of dual variables in continuum mechanics.” Continuum Mechanics 
and Thermodynamics 1 (1989), pp. 165-196. 

[92] D. L. Henann and L. Anand. “A large deformation theory for rate-dependent elastic—plastic materials with combined 
isotropic and kinematic hardening”. International Journal of Plasticity 25.10 (2009), pp. 1833-1878. 

[93] D. L. Henann and L. Anand. “Fracture of metallic glasses at notches: effects of notch-root radius and the ratio of the 
elastic shear modulus to the bulk modulus on toughness”. Acta Materialia 57.20 (2009), pp. 6057-6074. 

[94] D. L. Henann and L. Anand. “A large strain isotropic elasticity model based on molecular dynamics simulations of a 
metallic glass”. Journal of Elasticity 104.1-2 (2011), pp. 281-302. 


41 


[95] G. R. Hencky. Obituary of Gerhard R. Hencky, son of Heinrich Hencky. Published in San Francisco Chronicle, available 
at http://www.legacy.com/obituaries/sfgate/obituary.aspx?pid=169558738. 2014. 

[96] H. Hencky. “Uber den Spannungszustand in kreisrunden Platten mit verschwindender Biegungssteifigkeit”. Zeitschrift 
fur Mathematik und Physik 63 (1915), pp. 311-317. 

[97] H. Hencky. “Uber die Beziehungen der Philosophie des ,,Als Ob” zur mathematischen Naturbeschreibung”. Annalen 
der Philosophie 3 (1923). available at https://www.uni-due.de/imperia/md/content/mathematik/ag_neff/hencky_ 
als_ob.pdf pp. 236-245. 

[98] H. Hencky. “Uber einige statisch bestimmte Falle des Gleichgewichts in plastischen Korpern”. Zeitschrift fur Angewandte 
Mathematik und Mechanik 3.4 (1923), pp. 241-251. 

[99] H. Hencky. “Uber die Form des Elastizitatsgesetzes bei ideal elastischen Stoffen”. Zeitschrift fur technische Physik 9 
(1928). available at www.uni-due.de/imperia/md/content/mathematik/ag_neff/henckyl928.pdf, pp. 215—220. 

[100] H. Hencky. “Das Superpositionsgesetz eines endlich deformierten relaxationsfahigen elastischen Kontinuums und seine 
Bedeutung fur eine exakte Ableitung der Gleichungen fur die zahe Fliissigkeit in der Eulerschen Form”. Annalen der 
Physik 394.6 (1929). available at https://www.uni-due.de/imperia/md/content/mathematik/ag_neff/hencky_ 
superpositionl929.pdf pp. 617-630. 

[101] H. Hencky. “Welche Umstande bedingen die Verfestigung bei der bildsamen Verformung von festen isotropen Korpern?” 
Zeitschrift fur Physik 55 (1929). available at www.uni-due.de/imperia/md/content/mathematik/ag_neff/henckyl929. 
pdf, pp. 145-155. 

[102] H. Hencky. “The law of elasticity for isotropic and quasi-isotropic substances by finite deformations”. Journal of Rheology 
2.2 (1931). available at https://www.uni-due.de/imperia/md/content/mathematik/ag_neff/henckyjrheology31.pdf, 
pp. 169-176. 

[103] H. Hencky. “The elastic behavior of vulcanized rubber”. Rubber Chemistry and Technology 6.2 (1933). available at 
https://www.uni-due.de/imperia/md/content/mathematik/ag_neff/hencky_vulcanized_rubber.pdf pp. 217-224. 

[104] N. J. Higham. Matrix Nearness Problems and Applications. University of Manchester. Department of Mathematics, 
1988. 

[105] N. J. Higham. Functions of Matrices: Theory and Computation. Philadelphia, PA, USA: Society for Industrial and 
Applied Mathematics, 2008. 

[106] R. Hill. “On constitutive inequalities for simple materials - I”. Journal of the Mechanics and Physics of Solids 11 (1968), 
pp. 229-242. 

[107] R. Hill. “Constitutive inequalities for isotropic elastic solids under finite strain”. Proceedings of the Royal Society of 
London. A. Mathematical and Physical Sciences 314 (1970), pp. 457-472. 

[108] R. Hill. “Aspects of invariance in solid mechanics”. Advances in Applied Mechanics 18 (1978), pp. 1-75. 

[109] A. Hoger. “The stress conjugate to logarithmic strain.” International Journal of Solids and Structures 23 (1987), 
pp. 1645-1656. 

[110] R. Hooke. “Lectures de potentia restitutiva, or of Spring, Explaining the Power of Springy Bodies (1678)”. Early Science 
in Oxford, Volume VIII: The Cutler Lectures of Robert Hooke. Ed. by R. T. Gunther. Oxford University Press, 1931. 

[111] H. Hopf and W. Rinow. “Uber den Begriff der vollstandigen differentialgeometrischen Flache”. Commentarii Mathematici 
Helvetici 3.1 (1931), pp. 209-225. 

[112] T. Hughes and J. E. Marsden. “Some applications of geometry in continuum mechanics”. Reports on Mathematical 
Physics 12 (1977), pp. 35-44. 

[113] J. W. Hutchinson and K. W Neale. “Finite Strain JVDeformation Theory”. Proceedings of the IUTAM Symposium on 
Finite Elasticity (1982). Ed. by D. E. Carlson and R. T. Shield, pp. 237-247. 

[114] A. Imbert. Recherches theoriques et experiment ales sur V elasticity du caoutchouc, available at www.uni-due.de/imperia/ 
md/content/mathematik/ag_neff/imbert.rubber.pdf Lyon: Goyard, 1880. 

[115] D. F. Jones and L. R. G. Treloar. “The properties of rubber in pure homogeneous strain”. Journal of Physics D: Applied 
Physics 8.11 (1975), p. 1285. 

[116] J. Jost. Riemannian Geometry and Geometric Analysis (2nd ed.) Springer, 1998. 

[117] G. R. Kirchhoff. “Uber die Gleichungen des Gleichgewichtes eines elastischen Korpers bei nicht unendlich kleinen 
Verschiebungen seiner Theile”. Sitzungsberichte der Mathematisch-Naturwissenschaftlichen Classe der Kaiserlichen 
Akademie der Wissenschaften in Wien IX (1852). 

[118] J. Lankeit, P. Neff, and Y. Nakatsukasa. “The minimization of matrix logarithms: On a fundamental property of the 
unitary polar factor”. Linear Algebra and its Applications 449 (2014), pp. 28-42. DOI: 10.1016/j .laa.2014.02.012 

[119] M. Latorre and F. J. Montans. “On the interpretation of the logarithmic strain tensor in an arbitrary system of repre¬ 
sentation”. International Journal of Solids and Structures 51.7 (2014), pp. 1507-1515. 


42 


[120] M. Latorre and F. J. Montans. “Response to Fiala’s comments on ’On the interpretation of the logarithmic strain tensor 
in an arbitrary system of representation’”. International Journal of Solids and Structures 56-57 (2015), p. 292. 

[121] H. Le Dret and A. Raoult. “The quasiconvex envelope of the Saint Venant-Kirchhoff stored energy function”. Proceedings 
of the Royal Society of Edinburgh: Section A Mathematics 125.06 (1995), pp. 1179-1192. 

[122] T. Lehmann. “Anisotrope plastische Formanderungen”. Romanian Journal of Technical Sciences. Applied Mechanics 
17 (1972), pp. 1077-1086. 

[123] G. W. Leibniz. “Letter to Jacob Bernoulli, September 24, 1690”. Leibniz: Samtliche Schriften und Briefe, Reihe III: 
Mathematischer, naturwissenschaftlicher und technischer Briefwechsel. Ed. by H. HeB. Vol. 4. Berlin: Akademie Verlag, 
1995. 

[124] P. Ludwik. Elemente der technologischen Mechanik. available at https://www.uni-due.de/~hm0014/ag_neff/ludwik. 
pdf Berlin: J. Springer, 1909. 

[125] A. I. Lurie. Nonlinear Theory of Elasticity. Elsevier, 2012. 

[126] C. Man. “Hartig’s law and linear elasticity with initial stress.” Inverse Problems. 14 (1998), pp. 313-319. 

[127] J. E. Marsden and T. Hughes. Mathematical Foundations of Elasticity. Courier Dover Publications, 1994. 

[128] R. J. Martin and P. Neff. “The GL(n)-geodesic distance on SO(n)”. in preparation (2016). 

[129] R. J. Martin and P. Neff. “Minimal geodesics on GL(n) for left-invariant, right-0(n)-invariant Riemannian metrics”, to 

appear in The Journal of Geometric Mechanics (2014). available at arXiv:1409.7849. 

[130] R. J. Martin and P. Neff. “Some remarks on the monotonicity of primary matrix functions on the set of symmetric 
matrices”. Archive of Applied Mechanics 85.12 (2015), pp. 1761-1778. DOI: 10.1007/s00419-015-1017-4 

[131] L. C. Martins and P. Podio-Guidugli. “A variational approach to the polar decomposition theorem”. Rendiconti delle 
sedute dell’Accademia nazionale dei Lincei 66.6 (1979), pp. 487-493. 

[132] L. C. Martins and P. Podio-Guidugli. “An elementary proof of the polar decomposition theorem”. The American 
Mathematical Monthly 87 (1980), pp. 288-290. 

[133] G. Merrill. “Biographical memoir George Ferdinand Becker”. Memoirs National Academy of Science, XXI (1927). 

[134] A. Mielke. “Finite elastoplasticity, Lie groups and geodesics on SL(d)”. Geometry, Mechanics, and Dynamics - Volume 

in Honor of the 60th Birthday of J.E. Marsden. Ed. by P. Newton, P. Holmes, and A. Weinstein. Springer New York, 

2002, pp. 61-90. 

[135] P. W. Mitchell. “Hencky’s remarkable equation”. Australian Geomechanics 44.4 (2009), p. 41. 

[136] M. Moakher. “Means and averaging in the group of rotations”. SIAM Journal on Matrix Analysis and Applications 
24.1 (2002), pp. 1-16. 

[137] M. Moakher. “A differential geometric approach to the geometric mean of symmetric positive-definite matrices”. SIAM 
Journal on Matrix Analysis and Applications 26.3 (2005), pp. 735—747. 

[138] M. Moakher and A. N. Norris. “The closest elastic tensor of arbitrary symmetry to an elasticity tensor of lower symme¬ 
try”. Journal of Elasticity 85.3 (2006), pp. 215—263. 

[139] M. Moakher and M. Zeraf. “The Riemannian geometry of the space of positive-definite matrices and its application to 
the regularization of positive-definite matrix-valued data”. Journal of Mathematical Imaging and Vision 40.2 (2010), 
pp. 171-187. 

[140] G. Montella, S. Govindjee, and P. Neff. “The exponentiated Hencky strain energy in modelling tire derived material 
for moderately large deformations”. to appear in Journal of Engineering Materials and Technology, Transaction of the 
ASME (2015). available at arXiv: 1509.06541. 

[141] A. E. Moyer. “Robert Hooke’s ambiguous presentation of “Hooke’s law””. Isis 68.2 (1977), pp. 266-275. 

[142] J. G. Murphy. “Linear isotropic relations in finite hyperelasticity: some general results”. Journal of Elasticity 86.2 
(2007), pp. 139-154. 

[143] S. Muller and S. J. Spector. “An existence theory for nonlinear elasticity that allows for cavitation”. Archive for Rational 
Mechanics and Analysis 131.1 (1995), pp. 1—66. 

[144] P. Neff. Mathematische Analyse multiplikativer Viskoplastizitat. Ph.D. Thesis, Technische Universitat Darmstadt, avail¬ 
able at http://www.uni-due.de/~hm0014/Download_files/neffdiss.ps Aachen: Shaker Verlag, 2000. 

[145] P. Neff. A new support for using Hencky’s strain measure in finite elasticity - Seminar at the International Research 
Center for Mathematics and Mechanics of Complex Systems, Cisterna di Latina, Italy. Invitation by F. dell’Isola. 
http://memocs .univaq. it/?p=4184 Accessed: 2015-03-05. 2013. 

[146] P. Neff, B. Eidel, and R. J. Martin. “The axiomatic deduction of the quadratic Hencky strain energy by Heinrich Hencky 
(a new translation of Hencky’s original German articles)”. arXiv:lf02.f027 (2014). 


43 


[147] P. Neff, B. Eidel, F. Osterbrink, and R. Martin. “A Riemannian approach to strain measures in nonlinear elasticity”. 
Comptes Rendus Mecanique 342.4 (2014), pp. 254-257. 

[148] P. Neff, B. Eidel, F. Osterbrink, and R. J. Martin. “The Hencky strain energy || \ogU\\ 2 measures the geodesic distance of 
the deformation gradient to SO(n) in the canonical left-invariant Riemannian metric on GL(n)”. Proceedings in Applied 
Mathematics and Mechanics 13.1 (2013), pp. 369-370. 

[149] P. Neff, A. Fischle, and I. Munch. “Symmetric Cauchy-stresses do not imply symmetric Biot-strains in weak formulations 
of isotropic hyperelasticity with rotational degrees of freedom.” Acta Mechanica 197 (2008), pp. 19-30. 

[150] P. Neff, J. Lankeit, and A. Madeo. “On Grioli’s minimum property and its relation to Cauchy’s polar decomposition”. 
International Journal of Engineering Science 80.0 (2014), pp. 209-217. ISSN: 0020-7225. 

[151] P. Neff and I. Munch. “Curl bounds Grad on SO(3).” ESAIM: Control, Optimisation and Calculus of Variations 14.1 
(2008), pp. 148-159. 

[152] P. Neff, I. Munch, and R. J. Martin. “Rediscovering G. F. Becker’s early axiomatic deduction of a multiaxial non¬ 
linear stress-strain relation based on logarithmic strain”, to appear in Mathematics and Mechanics of Solids, doi: 
10.1177/1081286514542296 (2014). available at arXiv: 1403.4675. DOI: 10.1177/1081286514542296 

[153] P. Neff and I.-D. Ghiba. “The exponentiated Hencky-logarithmic strain energy. Part III: Coupling with idealized isotropic 
finite strain plasticity”. Continuum Mechanics and Thermodynamics 28.1 (2016), pp. 477-487. DOI: 10.1007/s00161- 
015-0449-y 

[154] P. Neff, I.-D. Ghiba, and J. Lankeit. “The exponentiated Hencky-logarithmic strain energy. Part I: Constitutive issues 
and rank-one convexity”. Journal of Elasticity 121.2 (2015), pp. 143-234. DOI: 10.1007/sl0659-015-9524-7 

[155] P. Neff, J. Lankeit, I.-D. Ghiba, R. J. Martin, and D. J. Steigmann. “The exponentiated Hencky-logarithmic strain 
energy. Part II: coercivity, planar polyconvexity and existence of minimizers”. Zeitschrift fur angewandte Mathematik 
und Physik 66.4 (2015), pp. 1671-1693. DOI: 10.1007/s00033-015-0495-0 

[156] P. Neff, Y. Nakatsukasa, and A. Fischle. “A logarithmic minimization property of the unitary polar factor in the 
spectral and Frobenius norms”. SIAM Journal on Matrix Analysis and Applications 35.3 (2014), pp. 1132-1154. DOI: 
10.1137/130909949 

[157] A. N. Norris. “The isotropic material closest to a given anisotropic material”. Journal of Mechanics of Materials and 
Structures 1.2 (2006), pp. 223—238. 

[158] A. N. Norris. “Eulerian conjugate stress and strain”. Journal of Mechanics of Materials and Structures 3.2 (2008), 
pp. 243-260. 

[159] A. N. Norris. “Higher derivatives and the inverse derivative of a tensor-valued function of a tensor”. Quarterly of Applied 
Mathematics 66 (2008), pp. 725-741. 

[160] R. W. Ogden. “Compressible isotropic elastic solids under finite strain - constitutive inequalities”. The Quarterly Journal 
of Mechanics and Applied Mathematics 23.4 (1970), pp. 457-468. 

[161] R. W. Ogden. “On stress rates in solid mechanics with application to elasticity theory”. Mathematical Proceedings of 
the Cambridge Philosophical Society 75 (1974), pp. 303-319. 

[162] R. W. Ogden. Non-Linear Elastic Deformations. 1. Mathematics and its Applications. Chichester: Ellis Horwood, 1983. 

[163] A. Ohara, N. Suda, and S. Amari. “Dualistic differential geometry of positive definite matrices and its applications to 
related problems”. Linear Algebra and its Applications 247 (1996), pp. 31-53. 

[164] W. A. Oldfather, C. A. Ellis, and D. M. Brown. “Leonhard Euler’s elastic curves”. Isis 20.1 (1933), pp. 72-160. 

[165] M. Ortiz, R. A. Radovitzky, and E. A. Repetto. “The computation of the exponential and logarithmic mappings and their 
first and second linearizations”. International Journal for Numerical Methods in Engineering 52.12 (2001), pp. 1431- 
1441. 

[166] X. Pennec. “Emerging Trends in Visual Computing”. Ed. by F. Nielsen. Berlin, Heidelberg: Springer-Verlag, 2009. 
Chap. Statistical computing on manifolds: from Riemannian geometry to computational anatomy, pp. 347-386. 

[167] X. Pennec, P. Fillard, and N. Ayache. “A Riemannian framework for tensor computing”. International Journal of 
Computer Vision 66.1 (2006), pp. 41-66. 

[168] J. P. Poirier and A. Tarantola. “A logarithmic equation of state”. Physics of the Earth and Planetary Interiors 109.1 
(1998), pp. 1-8. 

[169] S. D. Poisson. Memoire sur I’equilibre et mouvement des corps elastiques. Paris: L’Academie des sciences, 1829. 

[170] W. Pompe and P. Neff. “On the generalised sum of squared logarithms inequality”. Journal of Inequalities and Appli¬ 
cations 2015.1 (2015), pp. 1-17. DOI: 10.1186/sl3660-015-0623-6 

[171] L. Prandtl. “Elastisch bestimmte und elastisch unbestimmte Systeme”. Beitrage zur Technischen Mechanik und Tech- 
nischen Physik. Springer Berlin Heidelberg, 1924, pp. 52-61. DOI: 10.1007/978-3-642-51983-3_8 


44 


[172] W. D. Reinhardt and R. N. Dubey. “Eulerian strain-rate as a rate of logarithmic strain”. Mechanics Research Commu¬ 
nications 22.2 (1995), pp. 165-170. 

[173] W. D. Reinhardt and R. N. Dubey. “Application of objective rates in mechanical modeling of solids”. Journal of applied 
mechanics 63.3 (1996), pp. 692-698. 

[174] J. G. Resetnjak. “Liouville’s conformal mapping theorem under minimal regularity hypotheses”. Sibirskii Matematich- 
eskii Zhurnal 8 (1967), pp. 835-840. 

[175] H. Richter. “Das isotrope Elastizitatsgesetz”. Zeitschrift fur Angewandte Mathematik und Mechanik 28.7/8 (1948). 
available at https://www.uni-due.de/imperia/md/content/mathematik/ag_neff/richter_isotrop_log.pdf, pp. 205- 
209. 

[176] H. Richter. “Verzerrungstensor, Verzerrungsdeviator und Spannungstensor bei endlichen Formanderungen”. Zeitschrift 
fur Angewandte Mathematik und Mechanik 29.3 (1949). available at https://www.uni-due.de/imperia/md/content/ 
mathematik/ag_neff/richter_deviator_log.pdf pp. 65-75. ISSN: 1521-4001. 

[177] H. Richter. “Zum Logarithmus einer Matrix”. Archiv der Mathematik 2.5 (1949). available at https://www.uni-due.de/ 
imperia/md/content/mathematik/ag_neff/richter_log.pdf, pp. 360-363. ISSN: 0003-889X. DOI: 10.1007/BF02036865 

[178] H. Richter. “Zur Elastizitatstheorie endlicher Verformungen”. Mathematische Nachrichten 8.1 (1952), pp. 65-73. 

[179] R. S. Rivlin. “On the definition of strain”. Some recent developments in rheology : based on the 19J9 Bristol Conference 
of the British Rheologist’s Club. Ed. by V. G. W. Harrison. London: United Trade Press, 1950, pp. 125-129. 

[180] P. Rougee. “A new Lagrangian intrinsic approach to large deformations in continuous media”. European Journal of 
Mechanics - A/Solids 10.1 (1991), pp. 15-39. 

[181] P. Rougee. “The intrinsic Lagrangian metric and stress variables”. Finite Inelastic Deformations - Theory and Ap¬ 
plications. Ed. by D. Besdo and E. Stein. International Union of Theoretical and Applied Mechanics. Springer, 1992, 
pp. 217-226. 

[182] P. Rougee. Mecanique des Grandes Transformations. Vol. 25. Springer, 1997. 

[183] P. Rougee. “An intrinsic Lagrangian statement of constitutive laws in large strain”. Computers and Structures 84.17 
(2006), pp. 1125-1133. 

[184] O. Sander. “Geodesic finite elements of higher order”, to appear in IMA Journal of Numerical Analysis (2015). preprint 
available at http://www.igpm.rwth-aachen.de/Download/reports/pdf/IGPM356_k.pdf 

[185] C. Sansour. “On the dual variable of the logarithmic strain tensor, the dual variable of the Cauchy stress tensor, and 
related issues.” International Journal of Solids and Structures 38 (2001), pp. 9221-9232. 

[186] C. Sansour and H. Bednarczyk. “A study on rate-type constitutive equations and the existence of a free energy function”. 
Acta mechanica 100.3-4 (1993), pp. 205—221. 

[187] J Schroder, P Neff, and V Ebbing. “Anisotropic polyconvex energies on the basis of crystallographic motivated structural 
tensors”. Journal of the Mechanics and Physics of Solids 56.12 (2008), pp. 3486-3506. DOI: 10.1016/j .jmps .2008.08. 
008 

[188] J. Schroder, P. Neff, and D. Balzani. “A variational approach for materially stable anisotropic hyperelasticity”. Inter¬ 
national journal of solids and structures 42.15 (2005), pp. 4352-4371. DOI: 10.1016/j .ijsolstr.2004.11.021 

[189] T. Sendova and J. R. Walton. “On strong ellipticity for isotropic hyperelastic materials based upon logarithmic strain”. 
International Journal of Non-Linear Mechanics 40.2 (2005), pp. 195-212. 

[190] B. R. Seth. Generalized strain measure with applications to physical problems. Tech. rep. Defense Technical Information 
Center, 1961. 

[191] R. T. Shield. “Inverse deformation results in finite elasticity”. Zeitschrift fur angewandte Mathematik und Physik 18.4 
(1967), pp. 490-500. 

[192] F. Sidoroff. “Sur les restrictions a imposer a l’energie de deformation d’un materiau hyperelastique”. Comptes Rendus 
de I’Academie des Sciences 279 (1974), pp. 379-382. 

[193] M. Silhavy. The Mechanics and Thermomechanics of Continuous Media. Texts and Monographs in Physics. Springer, 
1997. 

[194] M. Silhavy. “Rank 1 convex hulls of isotropic functions in dimension 2 by 2”. Mathematica Bohemica 126.2 (2001), 
pp. 521-529. 

[195] E. B. Tadmor and R. E. Miller. Modeling Materials: Continuum, Atomistic and Multiscale Techniques. Cambridge 
University Press, 2011. 

[196] R. I. Tanner and E. Tanner. “Heinrich Hencky: a rheological pioneer”. Rheologica Acta 42.1-2 (2003), pp. 93-101. 

[197] A. Tarantola. Elements for Physics: Quantities, Qualities, and Intrinsic Theories. Heidelberg: Springer, 2006. 

[198] A. Tarantola. “Stress and strain in symmetric and asymmetric elasticity”, available at arXiv-.0907.1833 (2009). 


45 


[1991 C. Truesdell. “Mechanical foundations of elasticity and fluid dynamics”. Journal of Rational Mechanics and Analysis 1 
(1952), pp. 125-300. 

[200] C. Truesdell. “Hypo-elasticity”. Journal of Rational Mechanics and Analysis 4.1 (1955), pp. 83—131. 

[201] C. Truesdell. “The simplest rate theory of pure elasticity”. Communications on Pure and Applied Mathematics 8.1 
(1955), pp. 123-132. ISSN: 1097-0312. DOI: 10.1002/cpa.3160080109 

[202] C. Truesdell. “Das ungeloste Hauptproblem der endlichen Elastizitatstheorie”. Zeitschrift fur Angewandte Mathematik 
und Mechanik 36.3-4 (1956), pp. 97-103. ISSN: 1521-4001. DOI: 10.1002/zamm. 19560360304 

[203] C. Truesdell and W. Noll. “The Non-Linear Field Theories of Mechanics”. Handbuch der Physik. Ed. by S. Fliigge. 
Vol. Ill/3. Heidelberg: Springer, 1965. 

[204] C. Truesdell and R. Toupin. “The Classical Field Theories.” Handbuch der Physik. Ed. by S. Fliigge. Vol. III/l. Heidel¬ 
berg: Springer, 1960. 

[205] C. Vallee. “Lois de comportement elastique isotropes en grandes deformations”. International Journal of Engineering 
Science 16.7 (1978), pp. 451-457. 

[206] C. Vallee, D. Fortune, and C. Lerintiu. “On the dual variable of the Cauchy stress tensor in isotropic finite hyperelas¬ 
ticity”. Comptes Rendus Mecanique 336.11 (2008), pp. 851-855. 

[207] B. Vandereycken, P.-A. Absil, and S. Vandewalle. “A Riemannian geometry with complete geodesics for the set of 
positive semidefinite matrices of fixed rank”. IMA Journal of Numerical Analysis 33 (2013), pp. 481-514. 

[208] J. P. Wilber and J. C. Criscione. “The Baker-Ericksen inequalities for hyperelastic models using a novel set of invariants 
of Hencky strain”. International Journal of Solids and Structures 42.5 (2005), pp. 1547-1559. 

[209] H. Xiao. “Hencky strain and Hencky model: extending history and ongoing tradition”. Multidiscipline Modeling in 
Materials and Structures 1.1 (2005), pp. 1—52. 

[2101 H. Xiao, O. T. Bruhns, and A. Meyers. “Logarithmic strain, logarithmic spin and logarithmic rate”. Acta Mechanica 
124.1-4 (1997), pp. 89-105. 

[211] H. Xiao, O. T. Bruhns, and A. Meyers. “Existence and uniqueness of the integrable-exactly hypoelastic equation r = 
A(tr D)I + 2 pD and its significance to finite inelasticity”. Acta Mechanica 138.1 (1999), pp. 31-50. 

[212] H Xiao, O. Bruhns, and A Meyers. “Hypo-elasticity model based upon the logarithmic stress rate”. Journal of Elasticity 
47.1 (1997), pp. 51-68. 

[213] H. Xiao and L. S. Chen. “Hencky’s elasticity model and linear stress-strain relations in isotropic finite hyperelasticity”. 
Acta Mechanica 157.1-4 (2002), pp. 51-60. 

[214] H. Xiao and L.-S. Chen. “Hencky’s logarithmic strain and dual stress-strain and strain-stress relations in isotropic finite 
hyperelasticity”. International journal of solids and structures 40.6 (2003), pp. 1455—1463. 

[215] E. Zacur, M. Bossa, and S. Olmos. “Multivariate Tensor-Based Morphometry with a Right-Invariant Riemannian Dis¬ 
tance on GL + (n)”. Journal of Mathematical Imaging and Vision 50 (2014), pp. 19-31. 

[216] P. A. Zhilin, H. Altenbach, E. Ivanova, and A. Krivtsov. “Material strain tensor”. Generalized Continua as Models for 
Materials. Springer, 2013, pp. 321-331. 


46 


A Appendix 

A.l Notation 

• R is the set of real numbers, 

• R + = (0, oo) is the set of positive real numbers, 

• R" is the set of real column vectors of length n, 

• R raxm is the set of real n x m-matrices, 

• 1 is the identity tensor; 

• X T is the transpose of a matrix X £ R nxm , 

• tr(X) = ^"=1 is the trace of X £ R” x ", 

• CofX is the cofactor of X G R" xn , 

• {X, Y) = tr(A' T F) = Y^i,i =l i s the canonical inner product on R nx ", 

. ||x|| = s the Frobenius matrix norm on R nxra , 

• sym X = |(X + X T ) is the symmetric part of X £ R raxn , 

• skewX = |(X — X T ) is the skew-symmetric part of X £ R nxn , 

• dev„ X = X — t tr(X) • 1 is the n-dimensional deviator of X € R nxn , 

• (X, F) Atj/ic>K = p (devr, symX, dev„ symh) + /u c (skew X, skew Y) + ^ tr(X) tr(F) is the weighted inner prod¬ 
uct on R” x ™, 

• IIAU C ,« = is the weighted Frobenius norm on R" xn , 

• GL(n) = {A £ R" xn | det A 0} is the general linear group of all invertible A £ M” xra , 

• GL + ( n ) = {A £ R nx ™ | det A > 0} is the identity component of GL(n), 

• SL(n) = {A £ R nxn | det A = 1} is the special linear group of all A £ GL(n) with det A = 1, 

• O (n) is the orthogonal group of all Q £ R" xn with Q T Q = 1, 

• SO (n) is the special orthogonal group of all Q £ O(n) with detQ = 1, 

• Synr(n) is the set of symmetric, real n x n-matrices, i.e. S T = S for all S £ Sym(n), 

• Sym + (n) is the set of positive definite, symmetric, real n x n-matrices, i.e. x T Px > 0 for all P £ 
Sym +(n), O/iel”, 

• gl(n) = R raxn is the Lie algebra of all real n x n-matrices, 

• so(n) = {W £ R nxn | W T = —W} is the Lie algebra of skew symmetric, real n x n-matrices, 

• st(n) = {X £ R" x ” | tr(X) = 0} is the Lie algebra of trace free, real n x n-matrices, i.e. tr(X) = 0 for all 
X e sl(n), 

• Q C R" is the reference configuration of an elastic body, 
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\hp = Dip is the first derivative of a differentiable function ip\ Q C 
gradient, 


often called the deformation 


curl v denotes the curl of a vector valued function v: R 3 —> R 3 , 

Curlp denotes the curl of a matrix valued function p: R 3 —>• R 3x3 , taken row-wise, 

p\ kl — > R" is a continuously differentiable deformation with X7p(x) £ GL + (n) for all x £ SI, 

F = \7p £ GL + (n) is the deformation gradient, 

U = y/F T F £ Sym + (n) is the right Biot-stretch tensor, 

V = VFF T £ Sym + ( n ) is the left Biot-stretch tensor, 

B = FF t = V 2 is the Finger tensor, 

C = F T F = U 2 is the right Cauchy-Green deformation tensor, 

F = RU = VR is the polar decomposition of F with R = polar(F) £ SO(n), 

Eq = log U is the material Hencky strain tensor, 

Eq = log V is the spatial Hencky strain tensor, 

51 = DpW(F) is the first Piola-Kirchhoff stress corresponding to an elastic energy W = W(F ), 

5 2 = F~ l Si = 2 DqW(C) is the second Piola-Kirchhoff stress corresponding to an elastic energy W = 
W{C) (Doyle-Ericksen formula), 


T = Si F t = A og yW(log V) 
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p. 116] is the Kirchhoff stress tensor, 


det F 


t is the Cauchy stress tensor, 


rpBiot _ jj _ DuW(U) is the Biot stress tensor corresponding to an elastic energy W = W(U), 
L = FF~ 1 is the spatial velocity gradient, 

D = sym L is the rate of stretching or spatial strain rate tensor, 

W = skew L is the spatial continuum spin. 
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A.2 Linear stress-strain relations in nonlinear elasticity 

Many constitutive laws commonly used in applications are expressed in terms of linear relations between 
certain strains and stresses, including Hill’s family of generalized linear elasticity laws (cf. Section 4.2.1) of 
the form 

T r = 2 p, E r + A tr (E r ) ■ 1 (A.l) 


with work-conjugate pairs ( T r ,E r ) based on the Lagrangian strain measures given in (1.3). A widely known 
example of such a constitutive law is the hyperelastic Saint- Venant-Kirch.hoff model 

S 2 = 2 p E\ + A tr (Ei) 1 = p {C — 1) + — tr(C — 1) • 1 

for r = 1 and Tj = S 2 , where S 2 denotes the second Piola-Kirchhoff stress tensor. Similarly, a number of 
elasticity laws can be written in the form 

T r = 2 p E r + A tr(£V) ■ 1 

with a spatial strain tensor E r and a corresponding stress tensor T r . Examples include the Neo-Hooke type 
model 

a = 2 p E\ + A tr (E\) 1 = p (B — 1) + — tr (B — 1) • 1 
for r = 1, where Tj = a is the Cauchy stress tensor, the Almansi-Signorini model 

a = 2 fiE-i + A tr(T_i) 1 = p (1 — B _1 ) + — tr(l — B _1 ) • 1 
for r — — 1 and T_ 1 = a. as well as the hyperelastic Hencky model 

r = 2 [i log V + X tr(log V) ■ t 


16 


for r = 0 and Tq = t. A thorough comparison of these four constitutive laws can be found in 
Another example of a postulated linear stress-strain relation is the model 

T Biot = 2 p log U + A tr(log U) ■ 1, 

where T Blot denotes the Biot stress tensor , which measures the 11 stress per unit initial area before deformation’’’ 


28 . This constitutive relation was first given in an 1893 article by the geologist G. F. Becker 18 152 , who 


deduced it from a law of superposition in an approach similar to that of H. Hencky. The same constitutive law 
was considered by Carroll 38 as an example to emphasize the necessity of a hyperelastic formulation in order 
to ensure physical plausibility in the description of elastic behaviour. Note that of the constitutive relations 
listed in this section, only the Hencky model and the Saint-Venant-Kirchhoff model are indeed hyperelastic 
(cf. £3] Chapter 7.4]). 


A.3 Tensors and tangent spaces 

In the more general setting of differential geometry, the linear mappings F, U, C, V, B and R as well as various 
stresses at a single point x in an elastic body Q are defined as mappings between different tangent spaces: 
for a point x G fl and a deformation ip, we must then distinguish between the two tangent spaces T X Q and 
T v ( x )ip(Cl). The domains and codomains of various linear mappings are listed below and indicated in Figure 
[19} Note that we do not distinguish between tangent and cotangent vector spaces (cf. 63)). 

F,R: Tjfl —> T lp ( x ')ip( f2), 

U, C : TJl -> T x n , 

V, B : T^ x) ip(n) T v{x) ip(n). 
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x(t) e \7ip.x(t) £ T ¥>(a:( t )) ^(fi) 

C/, C : T x fl —)• T x 0 V, B : 7^(12) ->• 


-f, i? : —> T^ x ^<p(Q) 

(two-point tensors) 

Figure 19: Various linear mappings between the tangent spaces T x £l and 


The right Cauchy-Green tensor C = F T F, in particular, is often interpreted as a Riemannian metric on 
f2; Epstein 61 p. 113] explains that “the right Cauchy-Green tensor is precisely the pull-back of the spatial 
metric to the body manifold ”, cf. 127 . If 12 and <p(f2) are embedded in the Euclidean space R ra , this connection 
can immediately be seen: while the length of a curve x: [ 0 , 1 ] —► 12 is given by f n y/(x,x) d t, where (•, •) is 
the canonical inner product on R”, the length of the deformed curve ip o x is given by (cf. Figure 19 1 


[ v/(s(¥>°i), d t = [ yj (F(x) x, F(x) x) d t = [ yj (C(x) x, x) d t. 

Jo Jo Jo 


The quadratic form g x (v, v) = (C(x) v, v) at x £ 12 therefore measures the length of the deformed line element 
Fv at ip(x) £ ip(fl). Thus locally, 


dist Euc i idi¥ ,(n)(</?(a;), <p(y)) = dist geodi o(a;, y ), 

where distEuciid^n)^^)) viv)) = ll < / 3 ( ;r ) ~ ¥>(y)|| is the Euclidean distance between tp(x),tp(y) £ </j(I2) 
and dist g eod,n(a:, y) denotes the geodesic distance between x,y £ 12 with respect to the Riemannian met¬ 
ric g x (v,w) = (C(x)v, w). 

Moreover, this interpretation characterizes the Green-Lagrangian strain tensor E\ = 1(C—1) as a measure 
of change in length : the difference between the squared length of a line element v £ T x Sl in the reference 
configuration and the squared length of the deformed line element F(x) v £ T v r x \(p(fl) is given by 

\\F(x) x|| 2 - ||u|| 2 = (C(x) v, v) - {v , v) = ((C(x) -l)v,v) =2 (Ei(x) v, v ), 


where ||. || denotes the Euclidean norm on R". Note that for F(x) = 1 + Vu(x) with the displacement 
gradient Vu(x), the expression ||F(a;)u || 2 can be linearized to 

||.F(x) u || 2 = ||(1 + Vu(x)) i ;|| 2 = ((1 + Vu(x)) v, (1 + Vu(x)) v ) 

= (v, v) + 2 (Vu(x) v, v) + (Vzt(x) v, Vu(x) v) 

= ||u || 2 + 2 (symVu(a;) v, v) + ||Vu(x) r '|| 2 
= ||t ;|| 2 + 2 (sym Vu(x) v, v) + h.o.t., 
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where h.o.t. denotes higher order terms with respect to Xu(x). Thus 

||P(x) v\\ 2 - ||w|| 2 = 2 (e(a;) v, v) + h.o.t., 
where £ = symVu is the linear strain tensor. 


A.4 Additional computations 

Let Cof F = (det F) ■ F~ T denote the cofactor of F £ GL + (n). Then the geodesic distance of Cof F to SO (n) 
with respect to the Riemannian metric g introduced in (3.51 can be computed directly by applying Theorem 

GED 

4istg eod (Cof F , SO(n)) 

dev n log \J (Cof F) T Cof F\\ 


= T 
= T 
= T 
= T 
= T 

= T 


K 

2 L 


tr (log \j (Cof F) T Cof F^ 

tr (log \J (detF) 2 • F~ 1 F~ T ^ 


K 

2 L 


dev n log y (detF) 2 • F~ 1 F~ T \\ 2 + - 
|dev n log \/F- l F- T \\ 2 + | tr (log ((det F) ■ l) + log Vf- 1 F~ t ^ 
|dev n log(C- 1 )|| 2 + | [tr ((Indet F) • 1 + log^" 1 ))] 5 
| — dev n log{/|| 2 + ^ [n ■ (In det U) — tr(logP)] 2 

dev n log P|| 2 + ^ ^ [tr(log U)] 2 . 


A.5 The principal matrix logarithm on Sym + (n) and the matrix exponential 

The following lemma states some basic computational rules for the matrix exponential exp: R nxn —► GL + (n) 
and the principal matrix logarithm log: Sym + (n) —> Sym(n) involving the trace operator tr and the deviatoric 
part dev ra X — X — 1:1 ^ • 1 of a matrix X £ R nxn . 

Lemma A.l. Let X £ R. nxn , P £ Sym + (n) and c > 0. Then 

i) det (exp (A)) = e tr(A) , 

_ tr(X) 

ii) exp(dev n X)=e " • exp(X), 

in) log(c • 1) = ln(c) • log(l), 

iv) log ((det P)" 1 /” • P) = log P - ln(d f P) • 1 = dev n log P. 

Proof. Equality i) is well known (see e.g. [22] ). Equality Hi) follows directly from the fact that exp: Sym(n) —> 
Sym + (n) is bijective and that exp(ln(c) • 1) = e ln ( c ) -1 = c ■ 1. Since AB = BA implies exp(AP) = 
exp(A) exp(P), we find 

tr(X) trfX) tr(x) 

exp(dev„ X) = exp(X-• 1) = exp(X) • exp(-• 1) = exp(X) • e » • 1, 

n n 

showing ii). For iv), note that 

tr(logP) = ln(detP) =3> log P - ln(d f P) ■ t = dev n log P , 

and 

exp(dev„ log P) = e • exp(log P) 

= (e ln(det ~ Vn ■ P = (det P)~ l / n -P. 

according to ii). Then the injectivity of the matrix exponential on Sym(n) shows iv). ■ 
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A.6 A short biography of Heinrich Hencky 

Biographical information on Heinrich Hencky, as laid out in 1196, 33 95 


November 2, 1885: Hencky is born in Ansbach, Franken, Germany 
1904: Hencky finishes high school in Speyer 
1904-1908: Technische Hochsclrule Miinchen 
1909: Military service with the 3rd Pioneer Battalion in Miinchen 

1912- 1913: Ph.D studies at Technische Hochscliule Darmstadt 
1910-1912: Work on the Alsatian Railways 

1913- 1915: Work for a railway company in Kharkov, Ukraine 
1915-1918: Internment in Kharkov, Ukraine 

1919- 1920: Habilitation at Technische Hochschule Darmstadt 

1920- 1921: Technische Hochschule Dresden 
1922-1929: Technical University of Delft 
1930-1932: Massachusetts Institute of Technology (MIT) 

1933-1936: Potato farming in New Hampshire 

1936-1938: Academic work in the Soviet Union, first at Kharkov Chemical Technical Institute, then at 
the Mechanics Institute of Moscow University 

1938-1950: MAN (Maschinenfabrik Augsburg-Niirnberg) in Mainz 

July 6, 1951: Hencky dies in an avalanche at age 65 during mountain climbing 



Hencky at MIT, age 45 


Hencky received his diploma in civil engineering from TH Miinchen in 1908 and his Ph.D from TH Darm¬ 
stadt in 1913. The title of his thesis was “Uber den Spannungszustand in rechteckigen, ebenen Platten bei 
gleichmafiig verteilter und bei konzentrierter Belastung” (“On the stress state in rectangular flat plates under 
uniformly distributed and concentrated loading”). In 1915, the main results of his thesis were also published 
in the Zeitschrift fur angewandte Mathematik und Physik 96 . 

After working on plasticity theory and small-deformation elasticity, he began his work on finite elastic 
deformations in 1928. In 1929 he introduced the logarithmic strain ei og = log ( o/ighian’m'th ) a tensorial 
setting 199 and applied it to the description of the elastic behavior of vulcanized rubber 1 03] . 

Today, Hencky is mostly known for his contributions to plasticity theory: the article “Uber einige statisch 
bestimmte Falle des Gleichgewiclits in plastischen Korpern” [98] (“On statically determined cases of equilib¬ 


rium in plastic bodies”), published in 1923, is considered his most famous work 196 
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